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Abstract 



We develop a cosmological model where primordial inflation is driven by a 'solid', 
defined as a system of three derivatively coupled scalar fields obeying certain symme- 
tries and spontaneously breaking a certain subgroup of these. The symmetry breaking 
pattern differs drastically from that of standard inflationary models: time translations 
are unbroken. This prevents our model from fitting into the standard effective field the- 
ory description of adiabatic perturbations, with crucial consequences for the dynamics 
of cosmological perturbations. Most notably, non-gaussianities in the curvature per- 
turbations are unusually large, with /nl ~ and have a novel shape: peaked in 
the squeezed limit, with anisotropic dependence on how the limit is approached. Other 
unusual features include the absence of adiabatic fluctuation modes during inflation — 
which does not impair their presence and near scale-invariance after inflation — and a 
slightly blue tilt for the tensor modes. 

1 Introduction 

There is certainly no shortage of models for primordial inflation. We regret to inform the 
reader that — as anticipated in the Abstract — we are going to add our own to this list. How- 
ever, we feel that the inflationary model we introduce here presents conceptually novel fea- 
tures that make it stand out as a radical alternative to the standard inflationary scenario. 
The main reason for this is that our model does not conform to the standard symmetry 
breaking pattern of inflationary models, and this has far-reaching implications. 

In the usual cases, the matter fields ipm feature time-dependent cosmological background 
solutions ipm{t), which spontaneously break time translations. As a result, there is one fluc- 
tuation mode TT^x) that can be identified with the associated Goldstone excitation. Roughly 
speaking, it can be thought of as an in-sync perturbation of all the matter fields, of the form 



When coupling to gravity is taken into account, such a mode describes adiabatic perturba- 
tions. As usual for Goldstone bosons, the spontaneously broken symmetry puts completely 
general, non-trivial constraints on these perturbations' dynamics [1]. This property is at 
the basis of the model-indepedent approach that goes under the name of "effective field 
theory of inflation" [2], whose tenets are particularly compelling since they encompass — at 
first glance — all cosmological models: cosmology is about time-dependent, homogeneous, 
and isotropic field configurations. 
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However, as we will see, there are other possibilities. In our case, we will be dealing with 
matter fields featuring time-independent, x-dependent background solutions. Apparently, 
this contradicts two facts about infiationary cosmology: 

1. The universe is homogeneous and isotropic; 

2. In an expanding universe physical quantities depend on time and, more to the point, 
that one needs a physical 'clock' — a time-dependent observable — to tell the universe 
when to stop inflating. 

As for item 1: x-dependent solutions can be compatible with the homogeneity and isotropy 
we want for cosmological solutions and for the dynamics of their perturbations, provided 
extra symmetries are imposed. For instance, to get an FRW solution for the gravitational 
field, we need an homogeneous and isotropic background stress-energy tensor. This can 
arise from matter fields that are not homogeneous nor isotropic, provided there are internal 
symmetries acting on the fields that can reabsorb the variations one gets by performing 
translations and rotations. The simplest example is that of a scalar field with a vacuum 
expectation value 

{(!)) = ax. (1.2) 

Such a configuration breaks translations along x. However, if one postulates an internal 
shift symmetry — )• + a, then the configuration above is invariant under a combined 
spacial translation/internal shift transformation. As we will see, this is enough to make the 
stress-energy tensor and the action for small perturbations invariant under translations. To 
recover isotropy as well, one needs more fields, and more symmetries. For instance — in fact, 
this is the case that we will consider in this paper — one can use three scalar fields (j)\x) 
(/ = 1, 2, 3), with internal shift and rotational symmetries 

+ = const , (1.3) 

(P' ^ 0'j4>' , O'j e S0{3) , (1.4) 

so that the background configurations 

(0^) = ax^ (1.5) 

are invariant under combined spacial translation/internal shift transformations, and under 
combined spacial/internal rotations. As we will review in sect. |2| such a system has the 
same dynamics as those of the mechanical deformations of a solid — the phonons. In this 
sense, the cosmological model that we are putting forward corresponds to having a solid 
driving infiation. If this interpretation causes the reader discomfort — in particular, if having 
a solid that can be stretched by a factor of ~ without breaking sounds implausible — one 
should think of our model just as a certain scalar field theory. As we will see, the structure 
of such a theory is the most general one compatible with the postulated symmetries — and 
the impressive stretchability we need can also be motivated by an approximate symmetry — 
so that from an effective field theory standpoint, ours is a perfectly sensible infiationary 
model. From this viewpoint, the fact that the solids we are used to in everyday life behave 
quite differently — quantitively, not qualitatively — seems to be an accident: they lack the 
'stretchability symmetry'. 
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Figure 1: The shape of non-gaussianities for our model, according to the standard conventions and 
definitions of ref. 



As for apparent contradiction number 2: In our model the role of the physical clock will 
be played by the metric. More precisely, it will be played by (gauge invariant) observables, 
made up of our scalars and of the metric, like for instance the energy density or the pressure. 
These can depend on time even for purely space-dependent scalar backgrounds, because 
in the presence of a non-trivial stress-energy tensor, the metric will depend on time, in 
a standard FRW fashion. Doesn't this correspond to a spontaneous breakdown of time 
translations too? At some level it is a matter of definition, but we will argue in sect. |4]that 
the operationally useful answer is 'no', in the sense that there is no associated Goldstone 
boson, and that one cannot apply to our case the standard construction of the effective field 
theory of inflation as given in [2j. 

Formal considerations aside, our peculiar symmetry-breaking pattern has concrete phys- 
ical implications, with striking observational consequences. For instance, it predicts a three- 
point function for adiabatic perturbations with a 'shape' that is not encountered in any 
other model we are aware of. Without going into details here, we display it in fig. [T] for the 
benefit of the non-gaussianity aficionados: it diverges in the squeezed limit, but in a way 
that depends on the direction in which one approaches the limit, with a quadrupolar angular 
dependence. Its overall amplitude is also unusually large, corresponding to /nl ~ 

Before proceeding to spelling out our model in detail, we close this Introduction with 
two qualifications. The first is that we have been using (and will be using) a somewhat mis- 
leading, but fairly standard, language: when spontaneously broken symmetries are gauged, 
the associated would-be Goldsone bosons are not in the physical spectrum — rather, they 
are 'eaten' by the longitudinal polarizations of the gauge fields. In fact, there is a gauge, 
the so-called unitary gauge, in which the Goldstone fields are set to zero. In our case we 
will be dealing with spontaneously broken translations and rotations, and when gravity is 
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dynamical, these are gauged. In unitary gauge one can set the scalars to their vevs (1.5), and 
have the corresponding excitations show up in the metric. The Goldstone language is still 
useful though, in that it captures the correct high energy/short distance dynamics of these 
excitations. For massive gauge theories, this statement goes under the name of "the equiv- 
alence theorem". For cosmological models, it is just the statement that at sub-cosmological 
distances and time-scales, in first approximation one can neglect the mixing between matter 
perturbations and gravitational ones. We hope the Goldstone boson nomenclature will be 
more useful than misleading. 

The second qualification is that our model is not entirely new. A different formulation of 
essentially the same inflationary model was put forward and briefly analyzed by Gruzinov in 
[3]. Our emphasis here will be on the peculiar symmetry breaking pattern, on the effective 
field theory viewpoint, on a systematic analysis of cosmological perturbations — including 
their non-gaussian features — and in general on the conceptual and technical differences with 
more standard inflationary models. Cosmological solids have also been used as an exotic 
model for dark matter |5], and more recently as dynamical media with negative pressure but 
well-behaved excitations [H]. 

As a general guide for the reader, the bulk of our (long) paper can be schematically 
divided into three parts: 



sects. |2]-[4] introduce our model and discuss its field-theoretical and conceptual aspects, 
including why it cannot fit into the standard EFT of inflation; 



sects. |5]-[8] contain a technical analysis of cosmological perturbations, up to their three- 
point correlation function; 



sects. [9] and [To] address a number of conceptual subtleties in the analysis of cosmological 
perturbations, stemming from a very unusual feature of our model: the absence of 
adiabatic perturbations during inflation. 



A brief summary of our results is contained in sect. 11, along with a number of possible 



generalizations of our scenario. Much technical material is contained in the Appendix. 



2 Effective field theory for solids (and fluids) 

To begin with, let us review how one can describe the mechanical degrees of freedom of a 
solid in modern effective field theory terms. The first systematic approach to this question 
has probably been that of ref. 0. Here instead we will review the equivalent construction 
of [H H] (see also [lO]), whose notation we will follow. 

Consider a medium, filling space. Let's neglect for the moment gravitational effects and 
let's take the metric to be flat. The medium's configuration space can be parameterized 
by its volume elements' individual positions. We can attach a (three-dimensional) comoving 
label (j)^ (with / = 1, 2, 3) to each volume element, and follow the volume element's trajectory 
in physical space, as a function of its comoving label and of time, x(0^, t). We can also do the 
opposite, which for our purposes turns out to be more convenient: since at any given time the 
mapping between 0^ and x is invertible, one can parameterize the system via (j)\x,t), that 
is, by keeping track of the comoving labels of the volume elements going through physical 
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position X, as a function of time. This is a completely equivalent description, yet it has the 
advantage of implementing the spacetime symmetries in the usual way: (x, t) = are just 
the usual Minkowski coordinates, and the 0'^'s are three Poincare-scalars. Our problem then 
is reduced to constructing a relativistic low-energy effective field theory for three scalar fields 
in Minkoswki space obeying certain internal symmetries, which we will discuss in a moment. 
To avoid future confusion with standard cosmological nomenclature, in the following we will 
call the 0^'s 'internal' coordinates, and reserve 'comoving' for the standard FRW coordinates 
for when we introduce gravity. 

Now, apart from faithfully following the volume elements they are associated with, the 
internal coordinates are completely arbitrary. There is a coordinate system however that 
makes discussing internal symmetries — the symmetries that act on the (j)^ fields — particularly 



easy: it is the choice that reduces to (1.5) for a medium at rest, in equilibrium, in an 
homogeneous state at given external pressure. The parameter a measures the compression 
level, which can be dialed physically by changing the external pressure. Then, since we are 



supposed to be describing an homogeneous medium, but the expectation values (1.5) break 
translations, we are forced to postulate the existence of internal shift symmetries 

0^ 0^ + , = const , (2.1) 

so that there is an unbroken linear combination of spacial and internal translations that 
mimics the physical homogeneity we are trying to reproduce. In other words, thanks to the 



symmetry (2.1), the internal space of the medium is always homogeneous. Whether this 



homogeneity is realized in physical x-space depends on the medium's configuration, or state. 



For a state like (1.5), it is. This logic is, of course, the same as we were employing in the 



Introduction to motivate our cosmological model. 

We may want to impose rotational symmetries. The background configuration (1.5) 
breaks spacial rotations. For a crystalline solid, like e.g. one with a cubic lattice, a discrete 
subgroup of these should be preserved. This can be achieved by demanding an internal 
symmetry 

0^ ^ 0'j<p' (2.2) 
for a generic matrix O belonging to the desired subgroup of S'0(3), so that the background 



(1.5) is invariant under combined spacial and internal rotations. One can also impose full 
5*0(3) invariance, in which case one would be describing an isotropic solid with no preferred 
axes — a 'jelly'. For simplicity, in this paper we will consider this S'0(3)-invariant case only. 



We will come back to comment briefly on less symmetric cases in sect. 11 



So, we are led to consider the most general low-energy theory for three scalar fields 



obeying Poincare invariance and the internal symmetries (1.3) and (1.4). The shift (1.3) 
forces each field to appear with at least one derivative. Lorentz invariance forces these 
derivatives to contract among themselves. At lowest order in the derivative expansion, the 
only Lorentz-scalar, shift-invariant quantity is the matrix 

B^^ = d^cj)^ d^'ct)-^ . (2.3) 

We then have to construct 5*0(3) invariants out of this matrix. For a 3 x 3 matrix, there 
are only three independent ones, which we can take for instance to be the traces 

[5], [B^], (2.4) 
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where the brackets [. . . ] are shorthand for the trace of the matrix within. Alternatively, one 
could take the determinant, and two of the traces above. In the following, we will find it 
convenient to use one invariant — say [B] — to keep track of the overall 'size' of the matrix B, 
and to choose the other two such that they are insensitive to an overall rescaling of B, e.g. 

^-m. y=^^.. (2.5) 

The most general solid action therefore is 

S = Jd^xF{X,Y,Z) + ... (2.6) 

where F is a generic function that depends on the physical properties of the solid — e.g. its 
equation of state — and the dots stand for higher-derivative terms, which are negligible at 
low energies and momenta. 



As a self-consistency check, notice that the background configuration (1.5) solves the 
equations of motion for any value of a, which can then be thought of as being determined 
by the boundary conditions (the external pressure we alluded to above.) The eom for the 
action above are 

^^(^.<^'5^)=0- (2-7) 



For a linear configuration like (1.5), all terms in parentheses are constant, since they depend 
on the 0^s' first derivatives, which are constant. The eom are thus trivially obeyed. 

As a side note, a perfect fluid fits into the same general description. Only, it features 
infinitely many more internal symmetries, which express the physical fact that for a fluid 
one can move volume elements around in an adiabatic manner without paying any energy 
price — only their compression matters. By contrast, if one tries to displace a solid's volume 
elements, even without compressing or dilating them, one encounters stresses that want to 
bring the volume elements back to their rest positions. Mathematically, in our language this 
property of a fluid is expressed as the invariance of the dynamics under internal volume- 
preserving diffeomorphisms: 

0^^e'(0), det^ = l. (2.8) 

Of the 5*0(3) invariants above, only one particular combination survives: the determinant 
oiB, 

det B = l{[Bf - 3[B][B''] + 2[B^]) , (2.9) 

which is of course insensitive to multiplication of B by unit-determinant Jacobians. At lowest 
order in the derivative expansion, a fluid's action is thus 

^fluid = jd'^xF{detB) + ... (2.10) 

So, in this formalism, a fluid is just a very symmetric solid. 
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Back to the solid. The background configurations (1.5) spontaneously break some of our 



symmetries. There are associated Goldsone bosons, which are nothing but fiuctuations of 
the ^'^'s about such a background, 



(2.11) 



We get these fiuctuations' free action by expanding our action ( |2.6[ ) to second order in tt^. 
Using 

5" = {S^-^ + d^TT^ + d-^Tc^ + dyd^-K-^) , (2. 12) 

after integrating by parts and neglecting boundary terms we get 



S^S2 = jdS[- \FxX ■ ^2 + {iFxX + |(Fy + Fz)) {di 

+ (iFxxX' + §^{Fy + Fz)) (V-Tf) 



(2.13) 



where the subscripts stand for partial derivatives, which are to be evaluated at the back- 
ground values 



F ^ 1/3 , Z -> 1/9 . 



(2.14) 



These Goldstone excitations are the solid's phonons. For what follows it will be convenient 
to split the phonon field vf into a longitudinal part and a transverse one. 



VX7fL = 0, V-7fT = 0. 



(2.15) 



It is straightforward to extract the longitudinal and transverse propagation speeds from the 
phonon's action: 



2 2 FxxX^ s jFY + Fz) 

3 FxX 9 FxX 

in terms of which the quadratic action is simply 



2 {Fy + Fz) 

3 FxX 



S,= d'x{-\FxX) ^^-^^{d,Ti,f-{cl-cl){y 



TT 



(2.16) 



(2.17) 



If we expand eq. (2.6) to higher orders, we get the interactions among the phonons. 



The expansion is straightforward, but already at cubic order the result is quite messy, and 
not particularly illuminating. Below, we will display explicitly the cubic interactions in a 
particular limit, which yields some simplifications, and which is the physically relevant limit 
for an infiationary background. For the moment, it suffices to say that, by construction, the 
n-th order interaction terms will be schematically of the form (Stt)'", where the derivatives can 
be spacial or temporal, and the indices are contracted in all possible ways. The coefficients 
of these interactions terms — the coupling constants — will be given by suitable derivatives 
of F, evaluated on the background solution. Like for all derivatively-coupled theories, our 
interactions become strong in the UV, at some energy scale Astrong- For our theory to be 
predictive for cosmological observables, we will need this strong-coupling scale to be above 
the Hubble rate H, for the whole duration of infiation. 
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3 Inflation 



We can now allow for a cosmological spacetime metric and for dynamical gravity, which, 



operationally, is trivial: the index-contraction in (2.3) should be done via g'^'^ rather than 



T]^'^, and the measure in (2.6) should carry a \/—g- As usual, 'minimal coupling' corresponds 



to the most general coupling one can have between a matter system and gravity at lowest 
order in the derivative expansion. Then our solid's stress-energy tensor is 

1 bS dF 



(3.1) 



As to the scalar fields' background configuration, the in (1.5) should now be inter- 



preted as comoving FRW coordinates. The reason is that the FRW metric is invariant under 
translations and rotations acting on the comoving coordinates, and we want the l.h.s. and the 



r.h.s. in (1.5) to transform in the same way under the symmetries we are trying to preserve. 



We can also choose the normalization of the comoving coordinates to set the a parameter 
to one, so that from now on the background configuration is simply 



X 



(3.2) 



When computed on the background, the stress-energy tensor reduces to the standard 
Ti"^ = diag( - p,p,p,p), with 



P 



X 



(3.3) 



where the subscript X stands for partial derivative, and F and Fx are evaluated at the 
background values for our invariants: 



Y ^1/3, Z ^1/9 . 



(3.4) 



Notice that — by construction — X is the only invariant that depends on the scale factor; Y 
and Z were designed to be insensitive to an overall rescaling of B^"^ . This is the reason why 
only Fx appears in the pressure: for an FRW solution, the pressure is related to the response 
of the system to changing the scale factor, i.e., the volume. For a more general configuration. 



the stress-energy tensor (3.1) has a more complicated structure, which depends on Fy and 



Fz as well, which we report here for later use: 

2FyY 3FzZ 



T^u = g,uF-2d,c^'d,(P' (^(^Fx 



X 



X 



6'' + 



2FyB 

X2 



IJ 



IK uKJ 



+ 



3FzB'^B 



Now, in order to have near exponential inflation, we need 

H 



Via the Friedmann equations 



3M2/ 



H 



2M 



ip + p) 



(3.5) 



(3.6) 



(3.7) 



PI 



^We are defining the Planck scale as = (SttG) 
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and eq. (3.3), we can express e directly in terms of our Lagrangian F: 



1 Fx _ dlogF 

F dlogX ' ^ ^ 



where we used eq. (3.4) for the background value of X. We thus see that if we want our 
solid to drive near exponential inflation, we need a very weak X-dependence for F. Which 
is not surprising, since X is the only invariant that is sensitive to the volume of the universe: 
for inflation to happen, the solid's energy should not change much if we dilate the solid by 
~ e®°; this is only possible if the solid's dynamics do not depend much on X. 

This also suggests how to enforce the smallness of e via an approximate symmetry. Con- 
sider the scale transformation 

(f)^^X(j)\ A = const. (3.9) 
The matrix B^"^ changes by an overall factor, which affects X but not Y nor Z. Therefore, 



the smallness of Fx can be interpreted as an approximate invariance under (3.9): If F only 



depended on Y and Z, it would be exactly invariant under (3.9), and this would prevent 
quantum corrections from generating some X-dependence. If we start with a small Fx at 
tree level, the symmetry is only approximate, yet all further X-dependence generated at 
quantum level will be suppressed by the small symmetry-breaking couphng constant — Fx 
itself. 

Notice that in general, ordinary scale invariance — that acting on the spacetime coordi- 
nates as — )■ Ax^ — cannot be readily used as an approximate symmetry to enforce the 
smallness of symmetry-breaking parameters. For instance, one cannot solve the Higgs-mass 
hierarchy problem this way Moreover, it is generically anomalous, that is, broken by 
quantum effects, as clearly displayed by the running of coupling constants for interactions 
that are scale-invariant at tree level. Here instead, we are dealing with a purely internal 



symmetry — eq. (3.9) — which commutes with all spacetime symmetries. It has nothing to do 



with spacetime scale-invariance. It is on an equal footing with our other internal symmetries 



(1.3), (1.4), and, like those, is non-anomalous and can be used to constrain the structure of 



the Lagrangian. To avoid confusion, in the following we will refer to the symmetry (3.9) as 
'internal scale invariance'. 

Notice in passing that at lowest order in the derivative expansion, imposing internal scale 
invariance gives us full (internal) conformal invariance as a byproduct. The reason is the 
following. The invariants we are using in our lowest-order action are combinations of traces 
of the form [5"] . Under a generic internal diff 

0'^e'(0), (3.10) 

these traces transform as 

[5"] ^ [(J^J) ■ B ■ (J^J) ■ 5 • • • (J^J) ■ B] , (3.11) 

where J is the diff 's Jacobian matrix, and we have used the cyclicity of the trace. Now, by 
definition, conformal transformations are the subgroup of diffs that change the fiat metric 
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5^^ only by an overall scalar factor, that is, they have a (J^J) proportional to the iden- 
tity. So, under a general conformal transformation, our traces only change by an overall 
(0-dependent) factor. Our Y and Z combinations are insensitive to such a change, and as 
long as the Lagrangian only depends on those, it is invariant under all internal conformal 
transformations. It is interesting that even though our internal scale invariance has in princi- 
ple nothing to do with ordinary spacetime scale invariance, it shares with it the unavoidable 
company of special conformal transformations. In both cases, special conformal transfor- 
mations are not needed to close the symmetry group, yet they are respected by generic 
scale-invariant dynamics. In the spacetime symmetry case, there are fundamental reasons 
why that happens |12l [13]. In our case, it appears to be an accidental feature of the lowest 
order truncation of the derivative expansion. 

Back to physics. The smallness of e in our model does not come for free. To see this, 
notice that by dialing Fx we could go continuously through e = and end up with negative e, 
that is, positive H, which violates the null energy condition (NEC). But our EFT conforms to 
the hypotheses of the general theorem of [8] , which links NEC violations to such pathologies 
as ghost- or gradient-instabilities and superluminality. For small but positive e, we could be 
dangerously close to these pathologies]^ Our quadratic action for the phonons — eq. (2.13) — 
is quite explicit about this: the phonon's kinetic energy is suppressed by e, 

Jd^xe\F\-jf^ + ... (3.12) 

For very small e, this can in principle lead to two problems for our theory: 



Superluminality: the gradient energies in (2.13) are not explicitly suppressed by e, and 



as a consequence the propagation speeds (2.16) are formally of order 1/e, unless the 



numerators are also small. In an effective field theory like ours, with spontaneously 
broken Lorentz invariance, superluminal signal propagation is not necessarily an incon- 
sistency. However, it prevents the theory from admitting a standard Lorentz-invariant 
UV-completion [l^. We therefore feel that it should be avoided. 

Strong coupling: unless interactions are also suppressed by suitable powers of e — and it 
turns that they are not — a smaller kinetic energy means stronger interactions. This is 



obvious if one goes to canonical normalization for tt, by absorbing the prefactor in (3.12 ) 
into a redefined phonon field. Then inverse powers of e will show up in the interaction 
terms, thus signaling that the strong coupling scale of the theory is suppressed by some 
(positive) power of e. We have to make sure that this strong coupling scale is above 
H, for the whole duration of inflation. 



As for the former issue, notice first of all that the term proportional to Fxx in the 
expression for c| is forced to be close to —2/3. The reason is that not only do we need the 



^These considerations are irrelevant for standard slow-roll inflation, where the smallness of e is achieved 
via an approximately flat potential. One cannot play with the potential's slope and end up with positive 
H. Positive H would require flipping the sign of the inflaton's kinetic energy, which would of course entail 
ghost-instabilities. On the other hand, in our case the sign of H is the same as that of Fx, which is the 
Lagrangian parameter we are playing with to make e small. 



10 



'slow-roll' condition (3.6) for inflation to happen, we also need 

^ = ^ « 1 (3.13) 

for inflation to last many e-folds|^ This forces the second derivative Fxx also to be small. 
In particular, given eq. (3.8), and 



we get 

- -l + t-\r]. (3.15) 



So that, at lowest order in e and r], the propagation speeds (2.16) reduce to 

It is quite interesting that in this limit they depend on exactly the same [Fy + Fz) combi- 
nation; we see no obvious reason why this should be the case. As a result, the two speeds 
are not independent: they are related by^ 

4^|(l + ci). (3.18) 

We thus see that for both speeds to be sub-luminal, we need 

cl<\. (3.19) 

that is, positive {Fy + Fz) (recall that FxX = eF is negative, because F is). We do not want 
{Fy + Fz) to be too positive though — otherwise, we end up with negative squared speeds, 
that is, exponentially growing modes. We thus need {Fy + Fz) to flt into a small window, 

< (Fy + Fz) < |e|F| . (3.20) 

We were able to motivate the smallness of Fx, 

FxX = eF , (3.21) 

via an approximate symmetry. It would be desirable to do the same for this new combination 
of derivatives. One possibility is of course to say that all derivatives of F are small, that is 



..B^^^eF. (3.22) 



OF 

^In the computations that follow, we will assume all the slow-roll parameters to be of the same order of 
magnitude. 

^To all orders in e and 77, the exact relation is 

4 = !(l + ci-ie+|?/) . (3.17) 
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This formally corresponds to an approximate invariance under all internal diffs 

0' ^ , (3.23) 

that is, to the statement that the value of the Lagrangian does not depend much on the 
fields. Less formally, and more physically, it corresponds to saying that the bulk of the 
solid's energy density and pressure are dominated by a cosmological constant, which does 
not depend on the fields. Although this is of course a technically natural choice — having a 
large cosmological constant was never a problem — it is not particularly interesting. It would 
be more interesting to find a symmetry that allows large derivatives of F, 

Fy, Fzr^ F (3.24) 

but that — in the limit of exact symmetry — forces the combination Fy+Fz to vanish. Another 



possibility would be a symmetry that makes Fy + Fz saturate the upper bound in (3.20) 



that is, that makes c| vanish. This is not as unlikely as it sounds: for instance the perfect 



fluid action (2.10) — whose structure is protected by the volume-preserving diff symmetry- 



features vanishing propagation speed for the transverse phonons P, |9], as can be checked 



explicitly in (2.16), by using expression (2.9) for the determinant. We have not been able 



to find a symmetry that enforces the condition (3.20) while preserving (3.24). We have to 
take such a condition as an assumption, which might involve some fine tuning, but which is 
nonetheless consistent and necessary for the consistency of our inflationary solution. 

As for the strong coupling issue, we have to estimate the strong coupling scale Agtrong 
in our small e limit, and make sure that cosmological perturbations are weakly coupled at 



horizon crossing, that is, at frequencies of order H. Expanding the action (2.6) to all orders 
in TT we get interactions of the form 

fn ■ (dnr , (3.25) 

where fn is some typical derivative of F. In our case some combinations of derivatives are 
small, 

FxX ~ {Fy + Fz) ~ eF , (3.26) 

but we do not expect this to yield a substantial weakening of interactions. For instance, 
we will see below that in our approximation the coefficient weighing cubic interactions is 
Fy, which, as we argued above, can be as large as the background energy density, Fy ~ F. 
Assuming that Fy is a good estimate for the coefficients encountered in interaction terms, 
and assuming for the moment that both cl and ct are of order of the speed of light — so that 
there is no parametric difference between time- and space-derivatives — we can estimate very 
easily the strong coupling scale: We can go to canonical normalization for the kinetic term 

£2 ~ ei^ ■ (dnf {dn,Y , (3.27) 

so that the n-th order interaction becomes 

Cnr^Fy. {dnr ^ j^.idTT^ ■ (3-28) 
This is a dimension-2?T, interaction, weighed by a scale 

A„~(5^)*-- (3.29) 
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(recall that F and Fy have mass-dimension four.) If Fy is of the same order as F, this is 
simply 

Kn ~ F^'^ ■ , {Fy ~ F) , (3.30) 

which, for n > 3 and e ^ 1, is an increasing function of n. The lowest of all such scales — 
which defines the strong coupling scale of the theory — is thus that associated with n = 3: 

Astrong = A3 ~ F'/^e'/' , {Fy ~ F) . (3.31) 

If on the other hand Fy is much smaller that F, say of order eF, we get that all interactions 
are weighed by the same scale, which then defines strong coupling: 

Astrong ~ A„ ~ F^/^ ■ e^/^ , {Fy ~ eF) . (3.32) 

Either way, the strong coupling scale is a fractional power of e smaller than the scale asso- 
ciated with the solid's energy density. 

If the propagation speeds cl, ct are non-relativistic, the estimate of the strong coupling 
scale depends on the specific structure of the interaction terms, that is, on how many time- 
derivatives there are. In general, one may expect stronger interactions, i.e., lower strong- 
coupling scales for non-relativistic excitations (see e.g. [9] for a systematic analysis of this 



phenomenon in a different limit of our solid action.) Notice first of all that, because of (3.18) 
the transverse phonon speed ct is always relativistic. 



f < 4 < 1 . (3.33) 

So, our estimates above always work for the transverse phonons' self-interactions. For lon- 
gitudinal phonons with cl ^ 1, we can repeat the estimate using the cubic interaction, 
assuming this is still a good indicator of the strong coupling scale of the theory. Expanding 
( 2.6[ ) up to cubic order, and neglecting terms that are proportional to X-derivatives of F or 



to (Fy + Fz), we find 

^3 ~ Jd'x{- ^^Fy) ■ |l6 [dir]' - 36 [dir]' {[dir ■ dir^] + [{dirf]) 



+ 18[{d7ry] + 18[{d7Ty -dn']^ , (3.34) 

where {dTi)ij = diiTj is the matrix of spacial derivatives of tt, dir^ is its transpose, and the 
brackets stand for the trace. Notice in particular that there are no time-derivatives. For 
Cl <C 1, one can estimate the strong coupling scale via the following trick [9]. We can redefine 
the time variable as t — ?• t' /cl- Now in the kinetic energy term there is no hierarchy between 
time- and space-derivatives, 

52 ~ y d^'xeF ■ (tt^ - cl{Vrr)^) J d^x' eF cl ■ {d'rrf , (3.35) 

and we can apply the usual order-of-magnitude estimates as for relativistic theories. The 
cubic interaction becomes schematically 

53 ^ / d^x' — ■ {d'n)^ . (3.36) 

J 
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To get the lowest possible value for the strong-coupling scale — that is the most dangerous 
one — we can take Fy to be as large as possible, Fy ~ F. nGoing to canonical normalization 
and estimating the strong-coupling scale as above we get|j 

Pstrong ~ F'/\e'clY/' (3.38) 

This is the strong coupling momentum scale, or equivalently, the strong-coupling energy 
scale in units that are appropriate for our new t' variable. To convert to the original units 
of energy, we have to multiply by an extra c^: 

i^strong ~ F'/'ie'ciy/'' . (3.39) 

As we mentioned above, cosmological perturbation theory is under control only if 

£^strong > H . (3.40) 

Relating H and F via the Friedmann equation, ~ F/Mp^, we get a lower bound on the 
combination e ■ c\, 

e-ci» {H/M^,f\ (3.41) 

In principle our H can be several orders of magnitudes smaller than the Planck scale, in 
which case this bound is not particularly restrictive. Still, it is a nontrivial condition for the 
self-consistency of the perturbative computations we will perform. 

Once the strong-coupling danger is exorcised, large interactions are demoted (or pro- 
moted) from a problem to an exciting feature of the model: they imply huge non-gaussianities 
for our cosmological perturbations. As we will see in sect. [5| our non-gaussian signal is 
peaked on squeezed triangles, with the same size-dependence as the so-called local forms of 
non-gaussianity, but with a different angular dependence. The corresponding /nl parameter 
is of order l/(e ■ c|), which is a factor of 1/e bigger than what one finds — at the same value 
of cl — in single- field models with non-relativistic sound speed. 

A clarification is in order: we have been analyzing the viability of our model focusing on 
the phonons' dynamics, neglecting the background spacetime curvarture and the phonons' 
mixing with gravitational perturbations. Of course this is not entirely correct. As we men- 
tioned in the Introduction however, at energies much bigger than H, or equivalently, for 
time-scales much shorter than H~^, curvature and mixing have negligible effects, and in first 



approximation they can be neglected. Our conditions above, (3.20) and (3.41), should then 



be thought of as necessary and sufficient for our system to be well-behaved in the UV, at 



""There is a third possibihty where Fy 3> F. While this may not be a problem when analyzing the strong 
coupling scale of the theory, it is still very unnatural to have the background energy scale much smaller than 
the other scales in the theory. We already have one cosmological constant problem in cosmology, it would 
maybe be better if we tried not to introduce a second one. 

^If one were to repeat the same analysis for a more generic n-th order interaction, also weighed by Fy 



like eq. (3.34), and also involving spatial derivatives only, one would get 

Atrong'-i^'/'e^cf^ , (3.37) 
which, for small e and c/,, is minimized at n = 3. 
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very short distances and time scales. Our detailed analysis of cosmological perturbations in 
sect. [5] will confirm these results. 

We should also point out that although we will be using standard 'slow-roll' nomenclature 
for the conditions (3.6), (3.13) and for the associated perturbative expansion, nothing is 
'rolling' in our system, slowly or otherwise: our (p^ scalars are exactly constant in time. As 
usual however, the so-called slow-roll expansion really relies on the slowness of certain time- 
dependent observables like if, if, etc., which are well defined regardless of the presence of 
a rolling scalar. We will still use 'slow-roll' to refer to such a weak time-dependence, hoping 
that this will not cause confusion. As we emphasized, in our case the physical origin of this 
slowness is the near independence of the dynamics on X, which is, among our invariants, 
the only one that depends on time. Besides e and 77, in the following we will need one more 
slow-roll parameter, 

(3.42) 

which is small, because depends on time only via the Lagrangian's X-dependence. 

Finally, we should comment on why we are focusing on a solid rather than on a perfect 
fluid. First, since eventually we will be interested in quantum mechanical effects — as usual, 
quantum fluctuations will be the 'seed' for cosmological perturbations — we focus on a solid 
because we do not know yet how to consistently treat the perfect fluid effective theory as a 
quantum theory |9j. The problem has to do with the transverse excitations, which appear 
to be strongly coupled at all scales. Second, even forgetting about the transverse excitations 
and focusing on the longitudinal ones, we would not be able to keep those weakly coupled 
for many e- folds. As clear from (3.17), to have vanishing c\ (which is one of the defining 
features of a fluid) and small e, we need rj ~ —1. But, by definition, rj = e/{He), so that we 
need e = FxX/F to decrease by an order one factor over one Hubble time, i.e. to decrease 
like some order-one power of 1/a. F has to be nearly constant over many e-folds, which 
means that it is actually the numerator FxX that is tracking l/a{t). But it is precisely 
combinations like FxX that control the strong-coupling scale for longitudinal excitations in 
a fluid [9j, which means that we cannot have FxX decrease by exponentially large factors 
without making the system strongly coupled at frequencies of order H at some point during 
inflation. 



4 Physical 'clocks' and reheating 

Eventually we want our inflation to end and to be followed by a standard hot Big-Bang 
phase, that is, we want the universe to reheat and to become radiation dominated Q In our 
case, this process can be thought of as a phase-transition from a solid state to a relativistic 
fluid state. The advantage of our language in dealing with such a transition is that, as 



^It should be noted that [3] avoided going through a specific model of reheating by simply evolving the 
solid until ct 0, and demanding that the solid turn into a perfect fluid at that point. Given our analysis 



above, if we want to keep all the slow roll parameters small, as is reflected in (3.33) cannot become much 
smaller than 3/4. Via our EFT approach we thus see that the reheating model of [3] necessarily entails a 
breakdown of the slow-roll expansion before reheating, which is, of course, a consistent possibility. 
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we emphasized in sect. [2| it describes both sohds and fluids in terms of the same long- 



distance degrees of freedom, our scalars (p^ . Only, the fluid action (2.10) enjoys (many) more 
symmetries. So, regardless of the microscopic dynamics that are actually responsible for the 
phase transition, at long distances and time scales reheating corresponds to some sort of 
symmetry enhancement of our action. We will be more specific about this in a moment. 

In terms of our infrared degrees of freedom, what triggers reheating? In standard slow- 
roll inflation, it is the inflaton itself, when its time-dependent background field reaches a 
critical value. On the other hand, in the absence of perturbations, our 0^'s are exactly 
time- independent: (0^) = . However the metric is not, and there are gauge-invariant 
combinations like our 

X = g'^'' 8^(1)' 8,(1)' , (4.1) 



or the energy density and pressure in eq. (3.3 ), that do depend on time. Usually we are used to 
solids turning into liquids — that is, melting — when the temperature exceeds a critical value. 
But we can also envisage a solid that 'melts' at zero temperature, when one of the physical 
quantities above goes past a critical value. Helium offers an example of such a phenomenon: 
at zero temperature one can turn liquid helium into a solid by raising the pressure beyond 
~ 25 atm, and melt it back again by lowering the pressure below that value. In our case, 
we need this zero-temperature melting to be associated with a substantial release of latent 
heat, so that the fluid we end up with is (very) hot As far as we can tell, this does not 
violate any sacred principles of thermodynamics. 

Notice that as far as the background solution is concerned, to lowest order in the derivative 
expansion all choices for what observable triggers our solid's melting are physically equiv- 



alent: this is evident in our parameterization of the action (2.6), where all time-dependent 
observables depend on time only through X, or, equivalently, through a(t). However, when 
we include fluctuations, we break this equivalence. For instance, the three invariants [B], 
[B"^], and [B^] are independent combinations of the fields' derivatives. In the presence of 
fluctuations, the hypersurface defined by [B] reaching its critical reheating value is different 
from that defined by [B"^] or [B^] reaching their critical reheating values. As a result, some of 
our predictions for cosmological correlation functions might depend on the physical variable 
chosen to trigger reheating. Notice that after reheating, in the hot fluid phase, there is no 



ambiguity: the lowest-order action (2.10) only depends on one variable, the determinant 
of B'-' . It is thus natural, although not obviously necessary, to postulate that reheating is 
triggered by the value of det B. 

So, in terms of our action, our assumption is that for large det B the action has the 



general structure (2.6), whereas for detS below a critical value, the action has the more 



restricted form (2.10). In the space of our X, Y, Z invariants, this corresponds to dividing 



up the space into two regions, where the action has different symmetries: eqs. (1.3), (1.4) 



for the former, eq. (2.8) for the latter. Moreover, the slow-roll condition H ^ H during 



inflation is protected by another (approximate) symmetry of the solid phase, eq (3.9), which 
we want to be maximally violated in the post-reheating fluid phase, which has H ~ H^. 
Since renormalization is local in field space, the existence of different regions with different 
symmetries is protected by precisely those symmetries, and is thus a consistent and natural 



^In the fluid phase, which is described by the action ( |2.10 1, the temperature is given hy T — '^^ 



dVdot B 
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assumption |^ For illustrative purposes, consider for instance the following action: 

L-i'(A,r,Z)oc| 5)2/3 for det5<l, ^^'^^ 

where / is a generic function that evaluates to one for the background values F — )■ 1/3, Z — )■ 
1/9, and we suppressed an overall common factor, which defines the density at reheating. The 
'gluing' at det B = 1 can be smoothed at will. In the first regime, C describes a solid driving 
an inflationary phase with constant e (for fixed Y and Z, det B scales like — hence the 
e/3 power). In the second regime, £ describes an ultra-relativistic fluid, with p = p/3 oc T"^ 
[21 122] • The two regimes have different internal symmetries, as discussed above, but they 
share the same degrees of freedom. The classical evolution of the background solutions and 
of perturbations can then be followed smoothly through the transition region, as the long 
wavelength degrees of freedom are the same all along, and the equations of motion are regular 
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Notice that we have been implicitly assuming that reheating is instantaneous, that is, 
that our solid/fluid phase transition happens in a time interval that is much shorter than 
the Hubble scale, which is reasonable in principle, but not necessary. One can also consider 
much slower transitions, which in field space would correspond to replacing the sharp critical 
values we have been talking about for our observables, with much more continuous transition 
regions. All our physical considerations above apply unaltered. Since, as we will discuss, 
some of our predictions are potentially model-dependent, for what follows we need to assume 
a specific model for reheating. So, we will assume that reheating is fast, much faster than 
H, and that it is controlled by the value of det B. 

4.1 Why not the EFT of inflation? 

Before turning to a detailed analysis of cosmological perturbations, we close this section by 
discussing why our model does not conform to the standard EFT of inflation. As we just 
saw, we do have physical 'clocks', that is time-dependent background observables, so why 
can't we use the standard results for spontaneously broken time-translations? The reason is 
that these time-dependent observables depend on time only because the metric does. 

To see why this subtlety is important, consider first the dynamics of our system at very 
short distances — in the so-called decoupling limit — where the matter fluctuations decouple 
from the gravitational ones and the Goldstone boson language is appropriate. In first ap- 
proximation, this limit corresponds formally to setting G to zero. But without gravity, our 
background solution has no time-dependence whatsoever! All observables like density, pres- 
sure, etc. are now exactly constant in time, and only spacial translations and rotations are 
broken by the background configuration (0) = . As a result, the Goldstone bosons, whose 
existence and properties have to be assessed in the decouphng limit — because as recalled 



similar mechanism is at work in ghost inflation 23 , or in the EFT description of finite-temperature 
superfluids [24 . 
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As mentioned above, one should refrain from performing quantum computations in the fluid phase [9_ 



So, eq. (4.2) should not be thought of as a quantum effective theory in the second regime. Still, since at 



reheating all relevant modes are well outside the horizon and, thanks to the usual reasons, can be treated as 



classical, we only need (4.2 1 to be a consistent classical field theory, which it is. 
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in the Introduction, only in this hmit does it make sense to talk about them — are those 
associated with this spontaneous symmetry breaking pattern, not with time-translations. 

It is not surprising then, that once we re-introduce gravity, the dynamics of cosmological 
perturbations at all scales are quite different than for the EFT of inflation. This is manifest 
in the so-called unitary gauge, where one chooses the time-variable according to a physical 
clock. In the standard case, that clock would be the inflaton, and choosing the equal-time 
surfaces to be the equal-inflaton surfaces automatically sets to zero the inflaton perturbations 
and makes the metric the only fluctuating field. In our case, if we choose one of our time- 
dependent observables — say det B or [B] — to define unitary gauge, we are still left with 
matter perturbations, because the background time-dependence that we are using to set 
the gauge is carried by the metric, not by matter fields. We can either include the matter 
perturbations explicitly in the Lagrangian terms that we write down in this gauge, or we 
can set them to zero, by choosing the spacial coordinates now so that (p^ = . This is a 
complete gauge-fixing — all spacetime coordinates have been unambiguously defined — and is 
of course quite a different gauge choice than the standard unitary gauge. In particular, it is 
inequivalent to choosing the so-called ^-gauge for spacial diffs. Either way, the Lagrangian 
terms one would write down are quite different than for the standard EFT of inflation. We 
will go into the details of this new unitary gauge in sect. [B| 

The presence of matter fluctuations in the 'naive' unitary gauge cannot be taken as a 
sign that we are dealing with what would be called a multi-field model in the standard 
classification. First, because there is a gauge — our 'improved' unitary gauge — in which all 
matter fluctuations are set to zero. Second, because our spectrum of cosmological pertur- 



bations only includes one scalar mode, as clear from the Goldstone quadratic action (2.17). 
Furthemore, as we will see, this scalar mode is not adiabatic. In other words, in our system 
there are no adiabatic modes of fluctuation. This is yet another manifestation that we are 
dealing with a truly unconventional cosmological system. 

Notice that at the classical level, a subtle, isolated exception to all of the above is of- 
fered by a perfect fluid. On the one hand, in our language a perfect fluid is just a very 
symmetric solid. In particular, it features the same symmetry breaking pattern. On the 
other hand however, because of powerful conservation laws for vorticity, classically one can 
consistently set to zero the transverse excitations — the vortices — and be left with an EFT 
for the compressional modes only [8J. This admits a description in terms of a single scalar 
which spontaneously breaks tzme-translations — a P{X) theory — to which the standard EFT- 
of-inflation construction is applicable In particular, for a perfect fluid scalar cosmological 
perturbations are adiabatic. Once quantum effects are taken into account however, trans- 
verse excitations cannot be neglected any longer. In fact, for a fluid they are not particularly 
well behaved quantum mechanically [S], which is one of the reasons why we have been con- 
sidering a more generic solid rather than the special perfect-fluid case. 

As a technical aside, we should also emphasize that in a gauge where the matter fields 
are unperturbed, 0^ = x^, our B^"^ matrix reduces simply to g^"^ , and our Lagrangian thus 
becomes the sum of the Einstein-Hilbert action and of a particular function of g^"^ , that is, it 



^^Here we are using the standard notation of the community: P is a generic function and X here does not 
refer to the [B] we have been considering throughout the text but rather d^ipd^tp, where -0 is some scalar 
field whose background vev is (-0) = 
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reduces to a Lorentz- violating theory of massive gravity. Theories hke this have been studied 
in broad generahty in [25] . The reader famihar with the EFT of inflation might wonder why 
we are not writing down directly the action for the perturbations 5g^'^ in this gauge — the 
analog of unitary gauge in that case — , instead of going through the (apparently) unnecessary 
burden of writing an action for the full fields, solving for the background solution, and then 
expanding the action in small perturbations. The technical reason is that, unlike 5g^^ or 
6Kij for the EFT of inflation, our 5g^^ does not transform covariantly under the residual 
diffs, which are just time diffs for us. The reason is that g^"^ does, but its background value, 
j^^a^^^"' does not. It is then technically more convenient to write an action for the full g^'^ , 
which just amounts to writing an action for i?^"^, like we have done. 



5 Cosmological perturbations 



The three sections that follow contain a technical analysis of cosmological perturbations. 
Before skipping directly to sect. |8l the reader uninterested in the details of the derivations 



should be aware of our results: the scalar tilt (6.34), the tensor-to-scalar ratio (6.35), the ten- 



sor tilt (6.19), and the three-point function of scalar perturbations (7.15) (which is analyzed 
in some detail in sect. |8]). 

As the background stress tensor takes the usual homogeneous and isotropic form repre- 
sented by T^y = diag( — p,p, all the interesting repercussions of our pecuhar symmetry 
breaking pattern lie in the dynamics of perturbations around the slow roll background. In 
order to best isolate the dynamical degrees of freedom of the gravitational field it is most 
convenient to work in the ADM parametrization of the metric: 



(5.1) 



It is straightforward to check that the inverse metric g'^'^ is given by 
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where h^^ is the inverse spatial metric: h^^h^j = 6). For the background FRW metric 



9- 



ATS 



(5.2) 



0, and hij = a^it) 6ij. 



Following [18] we can write the action as 



S = jd^x NVh {iM^i ^ N-^^EijE'^ - E^)] + F{X, F, Z)} 



(5.3) 



where R^^"^ is the 3-dimensional Ricci scalar constructed out of hij and Eij = N Kij, with 
Kij denoting the extrinsic curvature of equal-time hypersurfaces. The constraint equations 



given by varying (5.3) with respect to A^ and A^* are: 



= iM^j _ j^-2^EijE'^ - E^)] + F{X, Y,Z) + N 
= iM^.V. [N-\E^ - 6]E)] + N ^^^fJ;^^ . 



dF{X,Y,Z) 
ON 



(5.4) 
(5.5) 
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The derivatives of F with respect to and A^-' can be calculated easily by noting that our 
B^^ (and hence X, F, Z) can be expressed in ADM variables as 

5" = - ^ (0' - Ar^'5,00 (0^ - N^d,(l)') + h^^^d.cl^'d^ct)' . (5.6) 

It is usually convenient to work in a gauge where scalar perturbations are removed from 
the matter fields and appear only in the metric, as gij = a^it) (1 + 2Q5ij, see e.g. [IH]. This 
possibility is not available to us: By using up the three spatial diffs, we can set the matter 
field perturbations to zero, (l)^ = , but then the spatial metric has an extra scalar mode, 
proportional to didjx, which we now cannot remove in the usual manner. However, we are 
still free to use time diffs, but these at best can set the scalar in front of 6ij, not x, to zero. 
A more useful choice is to use the time diff to set a physical "clock" — like those we discussed 
in the last section — to its unperturbed value. If this clock controls reheating, then reheating 
will happen at the same time for all observers in this gauge. We review this gauge choice, 
which we call 'unitary', in Appendix [B| 

For the moment we find it more convenient to work in spatially flat slicing gauge (SFSG) — 
defined in Appendix [B] — where we can write the fluctuations about the FRW background 
as 

(j)^ = x^ + 71^ , hij = a{tf exp(7ij) , N =l + 5N , (5.7) 
where is transverse and traceless, i.e. 

dilij = 7n = . (5.8) 

We can also further split the vr* and A^* fields in terms of their longitudinal scalar and 
transverse vector components. We therefore write: 

vr^ = ^=^L + 4 , and N' = ^=Nl + iV^ , (5.9) 

where 5i7r^ = diN^- = 0. From now on we will stop differentiating between internal I, J, . . . 
indices and spacial ones. The reason is that of the full original 5'0(3)spacctime x 

5*0 (3) internal Symmetry, only the diagonal combination is preserved by the background (j)^ = 
. TT* and A^* both transform as vectors under this unbroken 5*0(3), and therefore they 
carry the same kind of index. 

For our purposes here we are interested only in the leading non-gaussian behavior. Bar- 
ring accidental cancellations, this can be captured by keeping terms that are cubic in the 
fluctuations. In order to reproduce these terms it turns out to be necessary to only know A^ 



and A^' to first order in the fluctuations From now on, we find it easier to work in spacial 



Fourier space, with our convention defined for any field ^{x) by: 



^■^Tliis lucky fact is because the higher order terms in N and N'' will be multiplying the constraint 
equations. In particular: the third order term of N and TV' multiplies the zeroth order constraint equations, 
and the second order the first order constraint equations TF. If we where, however, to try and generate the 
fourth order terms we would need and iV* to second order. 
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For convenience, however, we will drop the twiddle as which field variable we intend will be 
obvious from the arguments. And so, solving the constraint equations (5.4) and (5.5) to first 
order in fiuct nations we have 

where the dot denotes a time-derivative. 

Now, plugging these solutions back into (5.3) will give us the correct action for the 
fiuctuations up to cubic order. For instance, the trilinear solid action after mixing with 
gravity is contained in Appendix D while the quadratic actions for the tensor, vector and 
scalar modes are contained in the next section. 

Now that we have the correct action for the perturbations in the presence of an infiating 
background we can compute correlation functions. In the end, we are interested in the post- 
reheating correlation functions of curvature perturbations, parameterized by either of the 
gauge invariant (at linear-order) combinations 

n = ^ + H5u, C = ^-H^- (5.14) 
2 2 p 



where we have followed the notation of [19] During our solid infiation phase, in spatially 
fiat slicing gauge these are given by 

where the non-local piece of TZ comes from solving the constraint equation for Nl. 

Two peculiarities concerning the behavior of these variables during solid infiation are 
worth mentioning at this point. First, TZ and C, do not coincide on super-horizon scales. 
Second, neither of them is conserved. These properties are in sharp contrast with what 
happens for adiabatic perturbations in standard cosmological models, and stem from the 
fact that during solid infiation, there are no adiabatic modes of fiuctuation! We will clarify 
why this is the case in sect. |9| 



^■^The general (i.e. before gauge fixing) perturbed metric (to the Unear-order) is parametrized by 

.g,j = a{tf {6,j{l + A) + d^djX + d,C, + d^Q + Aj) (5.15) 

with diCi — and diDij = Da — 0; furthermore the energy momentum tensor is decomposed into scalar, 
vector, and tensor modes as 

^Too = -pSgoQ + Sp (5.16) 

ST,o = pSgM ~{p + p){dM + 5uY) (5.17) 

5T,j = pSg.j + a^id^jSp + d.djSa + di6a' + djSa' + 5aJ,) . (5.18) 



21 



6 Two-point functions 



Upon plugging the expressions (5.11)-(5.13) back into the action, the quadratic action for 
tensor, vector, and scalar fluctuations reads: 
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translate into a (spacially) non-local structure in real space 



^\7rL-{H/H)7rL +Hci k' ^L 



(6.1) 
(6.2) 

(6.3) 
(6.4) 



(2) (2) 

Notice the quite nontrivial fc-dependence for Sj, and S)^ in Fourier space, which would 



6.1 Tensor perturbations 



Using (6.2) we can calculate the two-point function of the tensor perturbations. As usual, 
it is a simpler calculation than the scalar case and will serve as a warmup. We decompose 
the tensor modes into their polarizations 

^,,{k,t) = J2<jikh'{k,t), (6.5) 

s=± 

with efjefj* = 25^"^'. The transverse, traceless conditions on now simply become en = 
kiEij = 0. We further decompose each Y{k,t) as 

Y{k,t) = Yci{k,t) a%k)+YAk,tr a'\-k) . (6.6) 

— * — # 

where a^(fc)' and a^{k) are creation and annihilation operators obeying the usual commuta- 
tion relation 

[o^(^), a''\k')] = {27rf 6%k - k') 6''' , (6.7) 
and where the classical solution "^liik^t) obeys the equations of motion obtained by varying 

+ 2aH — 7e/ + [k^ + Ata'H''cl) 7d = . (6.8) 

In the above we have used conformal time r, where dr = dt/a, and the definition of the first 
slow-roll parameter ( |3.6[ ). Using the time-dependence of aH, e, and ct — which is worked 
out in Appendix [A] — we can express the e.o.m. for the tensor mode (6.8) up to first order in 
slow-roll parameters as 



T^^^i T^^i + + 5 

ar^ T dr \ r"^ 



Id - —^^Ici +\k' + 7d = . (6.9) 
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The subscript "c" denotes that the parameters c^, e are evaluated at some reference time r^, 
which is chosen to be the (conformal) time when the longest mode of observational relevance 
today exits the horizon, i.e. Tc is defined such that 

\cL,ckminrc\ ^ | CL.crc-f^today | = 1 (6.10) 

where the usual convention for flat spacetime atoday = 1 is understood. 
The most general solution to the above equation takes the form 

7.Kfc, r) = (-r)3/2+- [AHll\-kT) + BHi'Ji-kr)] , z.^ ^ | + 6, - |4,,e, , (6.11) 

where H^^''^'> are Hankel functions, and A and B are constants to be fixed by matching the 
appropriate initial conditions. 

At very early times the physical wavelength — k/a — is so small compared to the Hubble 
scale H that the curvature of spacetime cannot be perceived by such modes; it is therefore 
expected that the canonically normalized classical solution should match the free wave func- 
tion in the fiat-space vacuum, Note that the canonically normalized field 7can.(^) ''") 

is related to Yik, r) by 

Yik,r)^^l^. (6.12) 

Mpi a[T) 

Thus, the initial condition for 7c/(A;, r) is specified by 



Comparing this to the general solution given by (6.11) and using the asymptotic form (for 
large arguments) of Hankel functions, 

hm if«(x)-^./^e"-<"^+^)5, hm ^ /Ae--+<™+^)f (6.14) 

it is enforced that 

A = (_^j-ce^(-TV2+V4) + c>{e^), B = 0. (6.15) 

Where, once again. He = H{tc) and ec = e(rc). This result is valid up to first order in slow 
roll. 

It is interesting to note that even when these tensor modes are well outside the horizon, 
they are not conserved. A similar story applies to the gauge invariant curvature perturbations 
( and TZ defined by (5.19) and the vector perturbation (like vr^), as we will see in the following 
section. This is in opposition to the usual situation in most inflation models, and will be 
discussed in more detail in Section |9} In particular, by utilizing the asymptotic limit (for 
small argument) of the Hankel function 

\im H(^\x)-,H)^(-y , (6.16) 



23 



the mild time-dependence of the tensor mode in late time is given by: 



Jm ,7;,(J,r) = k-^''' (I)'"'"" f^^ + 0(,) ) . (6.17) 



where we have made use of relation (3.18). As we will see soon, the transverse vector modes 
and scalar modes in our model share this feature as well. 

And so finally, we are ready to obtain the two-point function for the tensor perturbations 
of the metric. In particular, we are interested in its late time behavior, when modes are well 
outside the horizon: 

{Y'ih,T)Y'ih,T)) = {2^f5\h + k2)5'''' |7d(^i,r)|' (6.18) 

The dependence on k and r is kept to first order in slow roll while the overall constant is to 
lowest order. 

The advantage of expressing time-dependent quantities in reference to a fixed fiducial 
time (tc), as opposed to the usual convention of using the time at horizon crossing, is that 
the time- and momentum-dependence are made manifest. We can simply read off the tilt of 
the spectrum to first order in slow roll from the above expression: 

nT - 1 ~ 2ci^,ee . (6.19) 

We can see that the two point function for tensor modes is blue shifted, which matches the 
result of [3], and which is a distinctive signature of our scenario, unreproducible by more 
conventional models of inflation. As to the spectrum's overall amplitude, it is the usual one: 
(77) ~ H^/Ml,. 

6.2 Scalar Perturbations 

We proceed by calculating the scalar two point function in a similar manner as above. As 



emphasized in section pi the scalar quantity of interest is the gauge invariant quantity C, 



which, in Fourier space, is related to the longitudinal Goldstones til simply by C = —kni/S 



(see eq. (5.19)). 

As before, let's decompose the scalar field of interest in terms of creation and annihilation 
operators: 

C(^,t) = Uk,t) bik) + Uk,tr h\-k) , (6.20) 
where the usual commutation relation is obeyed ^{k)^l>^ [k')] = (27ry^5^^\k — k'). 



^''We find C a more interesting quantity tlian TZ for the reason that given our assumption about reheating, 
C evolves continuously from inflation phase to post-inflation phase, while TZ does not. See Section [lO] for 
details. However, TZ does play a vital role as a simplifier in solving for the classical solution for scalar 
perturbation. 
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The classical equation of motion for C,c.i follows from varying the quadratic til action 



(6.4). The general solution to this equation is quite complicated, however there is a trick 
that makes its solution much easier. If we re-express the e.o.m. in terms of the other gauge 

(6.21) 



invariant parameter TZd (see eq. (5.19)) we have 
3ciCc/(fc,t) 



H 



n,i{k, t) + {3 + 7]{t) - 2e{t))n,i{k, t) . 



which, together with the definition of TZ 



n 



1 Cci + HeC,, 



d 



(6.22) 



He 1 + ky3amh ' 

forms a system of two first order equations of two variables. Eliminating we can generate 
a second order equation for TZ^i which takes the usual form similar to (6.8). Written with 
respect to conformal time, and up to first order in slow roll, it is given by: 

7^"^ + {2 + ri-2s) aH 7^;^ + [clk^ + (3e - 6s + Sc^e) a^H^] 7^d = , (6.23) 

where prime denotes a derivative w.r.t. conformal time, and s is the slow roll parameter 
defined by (3.42). Once again, using the conformal time dependence of aH, s, 1], e, and cl 
contained in Appendix |X] this equation can be solved in terms of Hankel functions. One 
finds that the general solution to first order in the slow roll parameters is given by 



7^c/(fc,r) 



CHl'Ji- cl{t) kT{l + s,)) + - CL{r) A;r(l + s,))] (6.24) 



where a = — ^(3 + 2ec + r/c — 2sc) and z^5 = ^(3 + 5sc — 2c\ + rjc). Notice that for this to 
be a solution, it is important to keep into account — to first order — the time-dependence of 
cl in the argument of the Hankel functions. 

Once again, in order to match the initial conditions we must canonically normalize ttj^. 
A quick glance at (6.4) reveals that the correct canonically normalized field is 

1/2 



T^rik,T) = ^'KL 



M^.a^k^ 

3 (l + 3 J me) 



-fcr— >oo 



2e MpiHa^ til ■ 



(6.25) 



With the usual normalization for the creation and annihilation operators we will recover the 
Minkowski vacuum for very early times by demanding that 



lim Cd{k,T) 

T— > — OO 

Or, equivalently that 

lim Tlci{k,T) 



3V2e MpxHa^ 



k 



-i{l + Sc)cL{T)kT 



AecL SMpiHa'^ 



a^H^ d 
k dr 



TT 



L,cl 

IT 



cl e 



-i(l + Sc)cL{T)kT 



Aek 



Mpia 



(6.26) 



(6.27) 



Matching the general solution given by (6.24) to the the initial condition (6.27) will set P = 
and 



C = —i. 



8e, Ml 



PI 



(6.28) 
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One can now use (6.21) to obtain the full expression for (ci{k,T), which is (as promised) a 
bit messy and not particularly instructive as for our computation we are only interested in 
C's late time limit. We will not bother to write it out here. 

Just like the tensor perturbations, neither TZ nor ( is conserved outside the horizon, 
though their temporal dependence is mild, i.e., suppressed by slow- roll parameters: 

lim 7^,Kfc,r) = (-) i-c^M^yl^^^-'^^/'-^^/' (- ^ ^\.^ ^ + 0(ei/2^ 

(6.29) 

4 2 / 

lim Cc/(^,r) = f-") (-c,,eA;r,)^i.--^-/2-''^/2 ( ^ , + 0{t'l'^ 



(6.30) 

Notice that at this order in slow-roll, on large scales ( and TZ are proportional to each other, 
with proportionality constant c|, which is in agreement with (6.21): 

7^~-4(r)C {kr^Q-). (6.31) 

Now finally, the two point function of ( for late times (when the modes are well outside 
the horizon) is given by 

(C(r,^i)C(r,fc2)> = i2nf5'ih + k2) \Ur,h)\' (6.32) 

-^^14° {2n)'6-'{k, + k2) X A . o2 (6.33) 

where, as before, we have kept the dependence on k and r to first order in slow roll while 
the pref actor is expressed only to lowest order in slow roll. 

Once again, since all the parameters are evaluated at the global time Tc as opposed to 
the time of horizon crossing for each mode, we can simply read off the tilt to first order in 
slow roll directly from the above expression. It is: 

ns~l^ 2eccl^^ - 5sc - Vc ■ (6.34) 

Notice the overall l/c| factor in front of the spectrum. In a more standard single-field model, 
this would be replaced by 1/cl (see e.g. [2]). For small cl, our extra powers of cl give us a 
very suppressed tensor-to-scalar ratio: 

r~eci. (6.35) 

It is crucial however to ascertain whether we should focus on the (( spectrum or the TZTZ 
one. After reheating, when the universe is dominated by a hot fiuid, they have to coincide, 
because of the usual reasons. But during infiation, because of (6.31 ), they differ by a factor of 
cj^ — precisely what suppresses our tensor-to-scalar ratio w.r.t the standard case. In sect. 10 
we argue that it is the (( spectrum that is continuous at reheating. 
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7 The three-point function 

We now compute the (CCC) three-point function. Like in single-field models with a small 
speed of sound, our three-point function will be enhanced by inverse powers of cl with respect 
to what one gets in standard slow-roll infiation, for essentially the same reason (see e.g. j2]). 
However in addition to this, we will find an extra 1/e enhancement, coming from the fact 
that the quadratic phonon action (3.12) is suppressed by e, whereas the cubic interactions 
(3.34) are not. 

In order to compute the correlation function at a specific time, we need to evolve it from 
a quantum state we know, that is the early-time fiat-space vacuum. Expanding the usual 
time-evolution operator and working to lowest order in perturbation theory we have the 
standard result, which is given schematically by: 



{arf) = -^ / dr' (^](-oo) | [C(r)^ H,^,{r')\ \ n{-oc)) . (7.1) 

For our purposes it is enough to calculate the three-point function to lowest order in slow-roll. 
As demonstrated in Appendix D, at this order it is enough to work with the phonon cubic 
action (3.34), which in our FRW curved background takes the form (neglecting boundary 
terms) : 

£3 = M|i a{tfH^ ^{^(dnf - IdTcdjTr'^dkTc' - pTidjU^djii^ + Idju'djTi^dkTt'^ . (7.2) 

Quite amazingly, this applies both in the decoupling limit {k ^ aHe^/"^) and in the opposite 
limit {k <^ aHt^/"^). And so, defining C^i = C(r, A;j), we have 

(CiC2C3> = ^Ml,^^-^ f {27rfS%p,+p2+ps)QipuP2,P3)y< (7.3) 
dr' a''(r')if^(r')([7rL(r, ki)nL{T, k2)nL{T, k^) , tx l{t\ Pi)ti l{t' , p2)'n l{t\ p-i)]) , 



where 

Q{PUP2,P3) = ^PlP2P3 - 7^ [Pl +P2 +P3 

81 27 V P2P3 pips P1P2 

^ 2 (Pl ■P2)(P2 ■P3)(P3 - Pl) .^^x 

3 P1P2P3 

which is totally symmetric under permutations of Pi,P2,P3- 

Writing vr^, in terms of creation and annihilation operators allows us to easily express the 
integral in terms of the classical solutions. To be precise, 

(CiCsCs) = -i^nfS'ih + k2 + h) X QMl,^^-^Q{h, h, h) Hr; -00) , (7.5) 
where the integral /(ri; T2) is defined as 

I{ri;r2) = Jin;T2) + J*{n;r2) (7.6) 
J{n;r2) ^ -znt{n,h)nt{n,k2)rft{ri,h) P dr' il^l4T^vr2'*(r', ^i)7r2'*(r', fc2)vrf *(r', 
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and we used that Q is an even function of the momenta. Just as in the previous section, we 
have used the dependence of H and a on conformal time with reference to Tc to lowest order 
in slow roll. 



Utilizing the classical solution of TZ{T,k) given by (6.24) we can immediately recover 



7r£ (r, k). This general solution, containing Hankel functions, and derivatives of Hankel func- 



tions is not particularly useful for attempting to perform the time-integral given by (7.6). 
We will first simplify 7rf (r, k), which will enable us to express the three-point function in an 
analytic form. 

First, note that at high momenta, above the de-mixing scale the quadratic Lagrangian 
for 7Ti takes the simpler form: 

^demix = jjdr [Ml.a^Hh) [|7r^|' - ciA;2|vri|'] , (7.7) 

which is — apart from the overall time-dependent pre-factor — the standard quadratic action 
for a scalar with generic propagation speed cl- Thus, one can solve the equations of motion 
generated by varying the above expression, fixing the exact form by requiring the flat-space 
vacuum in the infinite past (as we have done in the previous sections), without having to go 
through TZ. Of course, the two methods are equivalent, but this route makes the method of 
expansion clear. As for momenta that are deep in the full-mixing limit, where |c2,/cr| <C e^^^, 
TZ can be written in a simple form utilizing the asymptotic limit of the Hankel function 



(6.16), which upon insertion into (6.21) yields a particularly simple til = —<i(/k. In fact, 
even though in principle the two asymptotic expressions we get in this way should not share 
a common regime of applicability, in practice they are indistinguishable (at lowest order in 
slow-roll) in a wide range of momenta, e"*^*^^/"^-* < |cifcr| ^ 1. Therefore, as a very good 
approximation to the full solution, we can use the two asymptotic solutions and 'glue' them 
together anywhere inside this range. We find it more convenient to glue them at (c^/^r) ~ e. 
In summary, to lowest order in slow roll we can approximate the profile of 7r£' by: 

Bk{l + icLkT - |4A;V)e-^"^'=" , for |cifcr| > e(7.8a) 

Bk{ - CL,ckTf^-^''^'^{- CL,ckr,y^'^-'^^"~'^ + 0{e) , for |ciA:r| < 6(7.8b) 



<{r,k) 
where 



B, = -- ^ — (7 9) 



Our strategy now is to break up the integral (7.6) into separate regions where one of the 



functional forms described by (7.8a) and (7.8b) can used. The integral can then be done 
explicitly. 

7.1 Analytic Calculation of Integral 

To illustrate the point and make the flavor of the analysis transparent, let's look at an almost 
equilateral configuration of momenta. That is, assume that 

fci + A;2 + 4 = 0, fci ~ A;2 ~ A;3 ~ . (7.10) 
^'^This stems from the mildness of their time-dependence outside the horizon. 
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First, notice that given some reference time t^, we can split the time-integral as 

J[t;—oo) = zi; ^ zn— J[T^; —00) I' --Li; 



-2ec 



Then, choosing to be precisely the conformal time at which a mode of momentum k 
transitions from (7.8a) to (7.8b), —CLkr^ ~ e, we find that the real part of the second line 
vanishes at zeroth order in e, because all the vr^'s involved in the expression — inside and 
outside the integral — are real, and there is an overall i. The remaining piece is all that will 
contribute to the integral. And so we can write 



-5s/2-v/2-e^ 



J{t; -00) + J*{r; -00) = J] \Bkf {-CLhryl'+^-CLhr, 

i=l 

r dr ^m^^ n (1 - .c.^,r' - \clk]r-) e^-'^-' + c.c. 



1=1 



(7.12) 

where the scale invariant function ?7(fci, k2i k^) is given by 

+ 18(fc^fc2A;3 + kik^k^ + /^ifcaA;^) + 12(A;^A;^ + fc^fcf + klkl) 

+ 9(/i;fA;2 + 5 perms) + I2[klkl + 5 perms) + lS[klklk^ + 5 perms) | . 



In order to ensure convergence of the integral and project onto the right vacuum, the integral 
is actually computed over a slightly tilted contour, that is r' — )■ (1 — ie)T' + r*, with e — )■ 0+, 
and the limits of integration are from —00 to 0. Additionally, in the last step, the fact that 
1 ^ \cLkT^\ ~ e, and |r^,| > |r| has been used to collect only the leading order in slow roll 
contributions. 

A more careful analysis of the same flavor applies to more general triangle shapes, see 
Appendix [E| It turns out that the above expression is valid provided that the the triangle 
formed by the various momenta is not too squeezed, that is, provided 

fclong/^short > • (7.14) 

And so finally, putting everything together, we can express the full three-point function as 
(CiC2C3>(re) ^ (27r)353(^^ + ^2 + h) X (7.15) 
3 Fy Ht 1 (T,\^-l^ Q{k^MMu{k^,k2M) 



32 F M^^e^cl^ \rcJ kfklkl 
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where, we remind the reader, Q(ki,k2,k^) is given by (7.4) and U{ki, k2, k^) is given by 



(7.13) 



We will remark that the mild time dependence, (re/rc)^'^^'^, in the above expression can 
actually produce an order one correction to the overall magnitude of the three-point function. 
Indeed, assuming that inflation lasts for ~ 60 e-folds after the longest mode of today's 
relevance exits the horizon, we can see immediately that (Te/rc)'^^'^-' ~ £-60x0(6)^ which, as 
promised, depending on how small e is, can give an 0{1) correction . On the other hand. 



the mild momentum dependence {—CLkiTcYi-^ ^^/^ appearing in (7.12), is equal to one 



up to 0{e) corrections. We thus drop this piece from (7.15), in order to be consistent with 
the preceding computation of the integral. Our result (7.15) should be understood as the 
leading order contribution in slow roll. 

8 The size and shape of non-gaussianities 

It is useful to rewrite the three-point function above as an overall amplitude /nl times a 
shape that is a function of the momenta with order-one coefficients [4] . It is customary to do 
so at the level of correlators of the Newtonian potential $ during matter-domination, rather 
than of C- The relation outside the horizon is 

$ = |C. (8.1) 

Neglecting the tilt, one then writes the two- and three-point functions as 

mhmh)) = {2^f6\h + ^2) ■ ^ (8.2) 

(<l>(fci)$(fc2)<f (^3)) = {2Tif5\h + k2 + h) ■ fih, k2, ks) , (8.3) 
and normalizes / on 'equilateral' configurations with ki = k2 = k^ |16j . 



/(fc,A:,fc) = /NL-^. (8.4) 

This defines the parameter /nl in a model-independent fashion, in terms of observable 
quantities only. In particular, the observed value for the power-spectrum normalization 
is A$ ~ 2 X 10^^. Notice that, because of momentum conservation, the three momenta ki^2,z 
close into a triangle. As a result, the function / depends on the absolute values A;i_2,3 only, 
because a triangle is uniquely defined — up to overall rotations, which are a symmetry of 
/ — by specifying its sides. Notice also that scale-invariance forces / to have overall scaling 



dimension and we have used this fact in (8.4). Notice finally that the standard con- 
vention would be to call F the function that we call /. Unfortunately we have already been 
using F for our Lagrangian, so will stick to / for the function defined above. Hopefully this 
will not lead to confusion. 
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Applying these definitions to our case we find 



9 /tA**''" H 



'PI 

^2 . 



/(A:i,fc2,A:3) = --^^-^-J -A. ^pp3 (8-6) 

19415 /r,\-^^-^/^ Fy 1 Fy 1 

^^^ = -13122 -T-^'-^^'^-T-T^ ^'-'^ 

The /nl parameter gives us a measure of the absolute size of non-gaussianities. As we 
argued in sect. [3} Fy is essentially a free parameter, which can be as large as F, in which case 
our /nl is huge, of order l/(ec|). By comparison, single-field infiationary models with small 
sound speed — whose non-gaussianities are much larger than for standard slow-roll infiation — 
have, at the same value of the sound speed, an /nl which is a factor of e smaller than ours. 
Notice that in this case there is a potential tension for our model: the same combination e ■ 



appears in the scalar tilt, eq. (6.34). Of course one could have cancellations there, because 
of the other terms in the expression for the tilt. But assuming that these do not change the 
overall order of magnitude of the tilt, the tilt is small if non-gaussianities are large, and vice 
versa. Eventually, one should observe either. If on the other hand our Fy is of order eF, 
then this 1/e enhancement for non-gaussianities is gone, and our /nl becomes of order l/c|. 
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which is the same as for small sound-speed single-field models 

But the most interesting feature of our non-gaussian signal is probably its shape, that is, 
the dependence of / on the shape of the triangle made up by the momenta ^1,2,3- We plot it 
in fig. [T| following the standard conventions of [4J . In particular, it is clear from the plot that 
our three-point function is peaked on 'squeezed' triangles with <^ fci 2, but its behavior 
for those configurations depends strongly on the angle 9 between and the other momenta. 
Quantitatively, focusing on the ^^5^ structure in (8.6) and ignoring the prefactors from now 
on, we get 

40(l-3cos2e) , , 

/(/ci, k,, fca ^ 0) oc -- ^ ^3^3 ' . (8.8) 

where we used momentum-conservation to rewrite k2 as k2 — ki + k^ cos 9. Such an angular 
dependence is not there in any of the standard infiationary models we are aware of: at 
least for single-field models, the consistency relations force the behavior of the three-point 
function in the squeezed limit to be angle-independent — see e.g. [IHIE^- On the other hand, 
in our case the standard consistency relations are maximally violated, both at the level of 
angular dependence — as we just mentioned — and at the level of the overall prefactor: usually 
the squeezed limit of the three-point function is suppressed by the scalar tilt, which is of 
order e; here instead, there is no suppression like that, and in fact, as we argued above the 
overall prefactor can be as big as one over the tilt. It is easy to see why the consistency 



^^We remind the reader that cosmological perturbations can stih be nearly gaussian, even for huge values 
of /nl, as long as the combination /nl-\/ (C^) is much smaller than one at the relevant scales. For us, in 
the most strongly coupled case (Fy ^ F), such a combination is of order H/Mp\ ■ (ec|)^'^/^, which is much 



smaller than one if and only if the weak-coupling condition (3.41) is obeyed. As usual, perturbations are 



nearly gaussian if and only if they are weakly coupled at horizon crossing. 



31 



relations do not hold in our model: the standard argument of P^I2S] — that a long- wavelength 
background ( can be traded in for a rescaling of spatial coordinates — does not apply to our 
case, because in our model there is no gauge in which the curvature perturbation ( appears 
as a C ■ ^ij correction to the spatial metric. 

The fact that our non-gaussianities feature a novel shape, means that the numerical 
analyses of CMB data that have been carried out so far are quite suboptimal for our case. 
Following |3], this can be quantified by computing the overlap — or 'cosine' — between our 
shape and those used in the analyses. These cosines give a measure of how much one could 
improve by performing a dedicating analysis, which would entail using directly our shape in 
the estimator for the /nl parameter. Of all the standard shapes on the market, the /oca/ one 
is the only one that is peaked on squeezed triangles, with exactly the same scaling as ours, 
but with no dependence on the angle at which the squeezed limit is approached. As far as 



angular dependence goes, this corresponds to a monopole, whereas our squeezed limit (8.8) 
corresponds to a quadrupole. As a result, the overlap between our shape and the local one 
vanishes in the squeezed limit, and overall it is very small 

C0S(/, /local) ^ 2% . (8.9) 

Another popular shape is the equilateral one, which is peaked on configurations with ki = 
k2 = k^. Our shape is fairly suppressed for those configurations, and this results in a smallish 
overlap: 

C0S(/, /equil) ^ -39% . (8.10) 

Finally, given that our shape is fairly large for degenerate — or 'flattened' — triangles with 
ki = k2 + ks, one might think that it has a large overlap with the so-called orthogonal shape. 
However, we find a small cosine with the orthogonal template of [21] : 

cos(/,/orthog)^-32%. (8.11) 

All this suggests that (a) at the level of the three-point function, our model is very distin- 
guishable from more standard inflationary models, and that (b) to test our model, we should 
perform a dedicated CMB data analysis for our shape of non-gaussianities. 



8.1 A factorizable template 

Numerical analyses of three-point correlation functions in the CMB drastically simplify if 
one can rewrite the theoretical momentum-space three-point function in a factorizable form, 
that is, as a sum of a few terms that are just products of powers of the triangle sides fci, ^2, 
^3 [U]. If one is able to do so, the calculation time decreases by a factor of A^pixei, which 
is a huge improvement. Apart from isolated examples — most notably the so-called local 

^^According to the nomenclature of [4l, we are computing '3D' cosines, which are relevant for analyses of 
fully 3D data, like e.g. those of large scale structure surveys. For analyses of CMB data, which are projected 
onto a 2D sphere, one should refer to 2D cosines, which are however much more cumbersome to compute. 
Even though they are not entirely precise in this case, 3D cosines are still a good indicator of the overlaps 
between different shapes. We should also mention that for some of these cosines, the computations involve 
logarithmic divergences in the squeezed limit, which we cutoff at a value of 10'^ for the maximum ratio 
between momenta in the data set. 
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Figure 2: Our factorizable approximation to the three-point function. The flat triangular surface 
is the difference with our exact three-point function. As clear from the picture, this is maximized 
for equilateral configurations (x2 = X3 = 1), and is always quite small. 



form of non-gaussianities — factorization of a three-point function is typically possible only 
approximately. 

We present here a factorizable form that is a very good approximation to our three- 
point function. As for other forms of non-gaussianities, finding an adequate factorizable 
approximation is a matter of trial and error. The zeroth order desiderata are: (i) it should 
have an overall scale dependence; (ii) it should be totally symmetric in ki, k2, k^; 
(Hi) it should have the correct squeezed limit. Notice that our squeezed limit, eq. (8.8), is 
factorizable, since 

„ ki ■ ko , s 

cos^ = ^-^, 8.12 
kik2 

and the dot product can always be rewritten, using fci + ^2 + /cs = 0, in terms of absolute 
values, 

2fci . ^2 = fc^ - kl - kl . (8.13) 

So our first (and final) guess, is to take the squeezed limit behavior and rewrite it in a totally 
symmetric form: 

,2 -, 

(8.14) 



/factor(fcl, ^2, ^s) = " ^ 



kfkl-3(ki-h) 

h 2 permutations 



where the factor of 2 difference with eq. (8.8) takes into account that for /C3 ^ ki.,k2 we 
are getting the same contribution from the terms explicitly displayed and from one of the 
permutations. 

This ansatz obeys all the properties above, and performs surprising well in approximating 
our / away from the squeezed limit. Quantitatively, for flattened configurations {ki = 
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^2+^3) — where our three-point function has quite non-trivial features — the relative difference 
is very small, reaching a maximum of roughly 3% for 'folded' triangles (^2 = ^3 = ki/2). For 
equilateral configurations the difference is more substantial, but still quite small in absolute 
terms, given that our signal is small there. We plot this factorizable form in fig. |2| alongside 
its difference with the exact three-point function. The cosine between the two shapes is 

C0S(/, /factor) ^97%. (8.15) 



If this level of precision is not enough, one can improve (8.14) by adding to it the right 
admixture of the equilateral factorizable form of (TB], to make up for the small difference in 
the equilateral limit: 



/factor /factor(^l; ^2? ^s) + /equil(^l5 ^2; ^3) ^ ^ 



/cquil 

This way the cosine with / rises to roughly 99.5%. 



(8.16) 

k\=k2=ks 



9 Why Is Not Conserved? 

We saw in sect. [6] that already at linear level, during our solid inflation phase neither C, nor 
TZ is conserved on large scales. One might be tempted to attribute this to the presence of 
isocurvature modes in addition to adiabatic ones. However, in our model there is only one 
scalar perturbation — parameterized by vr^ in the gauge we have been using — and usually 
isocurvature modes are a luxury that only multi-field models can afford. To sharpen the 



paradox, our C, and TZ do not coincide on large scales — see eq. (6.31) — whereas usually they 
do, even for fluctuations that are not purely no n- adiabatic, that is, even when they are not 
conserved. In fact, Weinberg proved a no-go theorem stating that all FRW cosmological 
models — inflationary or not — feature two adiabatic modes of fluctuation, one of which has 
constant and identical C and TZ on large scales, while the other has C, = IZ = Q [T71 [TU] . 
This theorem is manifestly violated by our model. But, in Weinberg's own words, "no-go 
theorems have a way of relying on apparently technical assumptions that later turn out to 
have exceptions of great physical interest" [27j. We do not know whether our model will 
ultimately turn out to be of great physical interest, but it certainly offers an exception to an 
apparently technical assumption of ref. [T7j. Before showing this, let us explain in physical 
terms why our solid system cannot sustain adiabatic modes. 

By definition, an adiabatic mode is a perturbation that for very long wavelengths becomes 
locally unobservable, being indistinguishable from a slight shift in time of the background 
solution. An ordinary fluid offers a perfect example of this. Consider a long-wavelength 
sound wave in a fluid, for the moment in the absence of gravity. For an observer making 
measurements on scales much shorter than the wavelength, and working in the local rest 
frame — which is slightly different from the background one — , the only observables are the 
density and the pressure: neglecting gradients, a fluid is isotropic in its rest frame, and its 
stress-energy tensor is characterized by p and p only, which are related by the equation of 
state. Then, the only physical effect that is measurable on scales much shorter than the 
wavelength is the local compression (or dilation) the sound wave induces. When we include 
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gravity into the picture, essentially the same considerations apply for the perturbation, but 
now the time-evolution of the unperturbed FRW background already probes all possible 
compression levels for the fluid (within some range), that is, all possible values of p and p 
compatible with the equation of state. As a consequence, for wavelengths much longer than 
the Hubble scale, within any given Hubble patch a sound wave will be indistinguishable 
from a time-shift of the background solution, that is, it will become physically unobservable. 
Different Hubble patches will then evolve as separate identical FRW universes, and it can 
be shown that this translates into a conservation law for TZ [20] . 

For a solid, the situation is very different, already in the absence of gravity. A longitudinal 
phonon — which is the only scalar fluctuation at our disposal — does not correspond to a purely 
compressional deformation of the medium. Even for wavelengths that are much longer than 
the observation scale, the anisotropic stress and the compression associated with the phonon 
are of the same order of magnitude. This is evident from the form of the stress-energy tensor 



(3.5), which expanded to first order in vr yields schematically 

5T,j ~ (V ■ 7?) 5ij + {diTTj + djTT,) , (9.1) 

with similar coefficients in front of the two tensor structures — related to suitable derivatives 
of F w.r.t. X, Y, Z — whose precise values will not concern us here. For a longitudinal phonon 
of momentum q, the anisotropic stress is proportional to qiqj, and is of the same order 
of magnitude as the change in pressure. A local observer can detect this anisotropy if 
he or she can detect the change in pressure. In other words, unlike a fluid, a solid with 
small longitudinal deformations is not locally isotropic. Once we include gravity, these 
scalar fluctuations will be locally distinguishable from the background solution, even for 
super-horizon wavelengths, since the background is isotropic. One could probably apply a 
'separate universe '-like argument that includes anisotropic homogeneous backgrounds, to get 



information about the time-evolution of these modes, even at non-linear order ^ We leave 
this for future work. For the moment, we just notice that our scalar modes do not conform 
to the standard characterization of adiabatic perturbations. 

We can now go back to the no-go theorem of ref . [I7j , and see which assumptions we are 
violating. The theorem is based on the following ingenious idea. Since an adiabatic mode is, 
by definition, unobservable once the wavelength is very long, at zero momentum it should 
reduce to a gauge transformation of the FRW background. Newtonian gauge is a complete 
gauge-fixing at finite momentum, but it has a residual gauge freedom at zero momentum. By 
exploiting this gauge freedom one can construct zero-momentum solutions of the linearized 
perturbation equations. Most of these are pure-gauge, unphysical solutions. To be physical, 
they have to be the zero-momentum limit of finite momentum solutions, which are physical 
because there is no residual gauge- freedom at finite momentum. For this to be the case, one 
needs the zero-momentum solution to obey the finite-momentum version of the (ij) and (Oz) 
linearized Einstein equations, 

$ = ^ - (SttG) (5(T H6u = ^ + , (9.2) 

which singles out only two independent modes among the zero-momentum solutions. Here 



^^We thank Matias Zaldarriaga for this remark. 
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5a and 5u are the scalar anisotropic stress and the velocity potential of eqs. (5.17), (5.18), 

5Tij D a^{t) didjSa , SToi D -(p + p) d,5u . (9.3) 
One has to impose further that all equations of motion — for gravity and for the matter 



fields — are regular for g*— 0, including eq. (9.2). More precisely, if one rewrites all linearized 
equations of motion in first-order form. 



2/a(g, t) + Cab{q, t) ■ yb{q, t) = 



(9.4) 



where the i/aS include the fields and their velocities, with constraints for the initial conditions 
of the form 

Cbiq) ■ ybiq,to) = , (9.5) 

then one has to demand that all Cab and Cb coefficients be regular for q 0. This is the 
only technical assumption of the theorem If it is obeyed, then the two zero-momentum 
gauge modes can be promoted to physical, finite-momentum solutions, which are adiabatic 
by construction, and one of which turns out to have constant ( and TZ. 

From our physical argument above, we see that the culprit for us is the anisotropic stress, 
which does not become negligible at long wavelengths. Indeed, comparing (9.3) with (9.1) 
we get schematically 

g2 



6a 



TC 



(9.7) 



(we are neglecting factors of a(t), and corrections to (9.1) involving the metric, which does 
not change our conclusion.) Once plugged into eq. (9.2), this gives us an equation of motion 
that is not regular for g — )■ 0, thus violating Weinberg's technical assumption. We get a 
similar singularity in the second equation of (9.2), since from Toj ~ tt* we get a velocity 
potential 



Su 



q_ 

q2 



IT 



{9.t 



Notice that we cannot reabsorb the annoying q ^ factors into a new vr field, thus making 6a 



and 6u regular for g — )■ 0, because one of the equations (9.4) is the equation of motion for n 



itself, which is local in real space, thus analytic in q in Fourier space. If we were to divide it 
by g^, to write an evolution equation for the new vf/g^ field, we would introduce singularities 
there. 



^^There is also an implicit assumption — used to write down the zero-momentum pure gauge solutions — 
that all background matter fields only depend on time, which is not obeyed in our case. However, this is 
easily fixed by performing the correct gauge transformation, which, following the notation of ref. |17) , in our 
case yields the pure gauge solution 



7?= A 



Sp = —pe , Sp = —pe 



(9.6) 
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10 Matching of correlation functions at reheating, and 
post-inflationary evolution 

Since our scalar perturbations are not adiabatic, our predictions for post-inflationary cor- 
relation functions on large scales can in principle be affected by local physical processes 
happening at reheating. It would be interesting to investigate how much model-independent 
information about time-evolution across the reheating phase can be obtained by applying an 
anisotropic parallel-universe argument, as mentioned in the last section. For the time being, 
we adopt what appears to be a reasonable model for reheating. As motivated in sect. |4| we 
postulate that inflation ends when B = detB^"^ reaches some critical value B^,, after which 
the matter content of the universe turns into a perfect fluid, which can be described by a 
low-energy effective Lagrangian involving the same set of dynamical d.o.f.'s as our solid, but 
with different (and more restrictive) symmetry requirements. Furthermore, we assume that 
the transition (a.k.a. reheating) from solid to fluid occurs "smoothly" (in a sense that will be 
made explicit) and instantaneously (the transition time being much shorter than the Hubble 
time). 

Since it is the scalar quantity B that plays the role of the "clock" in our model, it is 
easiest to work in unitary gauge (UG), in which constant time slices correspond to surfaces 
of uniform B (the properties of this gauge are worked out in Appendix |B]). Explicitly, our 
model for reheating can be captured by the following statements: 

• Inflation ends at te, with a(te)~^ = Be, where t denotes the time in UG. 

• The matter content in the post inflationary era takes the form of a perfect fluid, which 
can be described by the effective action |9] 

5fluid = J d^x^F{B) (10.1) 

where B = detg^"^ in UG. The change from one equation of state to the other is 
effectively instantaneous. For an ultra-relativistic fluid with p = ^p, one has 

F{B) oc 5^/^ . (10.2) 

• Energy transfer from the solid phase to the fluid phase during this short reheating 
period is efficient, and the normalization of F is restricted in such a way that energy 
conservation is respected, i.e. 

Psoiid = -F{X,Y,Z)l^ = -F{a{te)-^) =pfl„id . (10.3) 

However, generally PsoUd 7^ Pfluid, since the equation of state has been changed. Conse- 
quently, even though the Hubble parameter H remains continuous, H does not. 

• Smoothness: the dynamical d.o.f.'s in unitary gauge — A (or x), Cj, Dij — as well as 



their first derivatives are continuous across the t = te surface ^ The second deriva- 
tives will exhibit discontinuities since the equations of motions are altered due to the 
instantaneous change in the equation of state. 

See Appendix U3\ for the definition of these fields in terms of fluctuations of the metric. 
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An immediate consequence following the smoothness requirement is that ( transits continu- 
ously from solid phase to post-inflationary fluid phase, while TZ does discontinuously. Indeed 
notice that in UG, ( and TZ are given by 

C = -, 7^ = -. ^ 10.4 

2' 2Hl-k^/?,a?H 

thus the discontinuity of TZ stems from that of H. 

Given that the transition from the solid phase during inflation to the perfect fluid phase 
during the post-inflationary era occurs effectively instantaneously, we can compute various 
correlators (of A, Ci and Dij) at and use them as the initial conditions of the post- 
inflationary evolution. There are two subtleties: 

1. How do we relate the correlators of quantities in UG to those in SFSG, which have 
been computed in Section |6] and Section [7} 

2. Given that the super-horizon modes of A and Dij are not adiabatic during inflation, 
they start the post-inflationary evolution with a non- vanishing first derivative in time. 
Although eventually these super-horizon modes become constant (time-independent) 
before reentering the horizon, a natural question to ask is how much the eventual 
constants could differ from the initial conditions these modes start with. 

In order to address the first issue, let's compute the scalar two- and three-point corre- 
lators: {Aiitf.) A2{te)) and (y4i(te) ^2(^e) ^3(^e)), whcrc the is to remind ourselves that 
we are using time in UG, and Aiit) is shorthand for A{ki, i). Using the transformation rule 
from SFSG to UG, we have that 

and we can write 

A{x) = + + . . . , 

where A'^'^\x) = f5i7r*(x) = 2C, A^^) ~ dirdir , etc. 

It follows immediately that, schematically, 

(A^) = (A^^M^i)) + 2(A«A(2)) + ... ~ 4(C') + {{dirf) + ... . (10.5) 

The second term on the right can be neglected, since it is of higher order in the perturbative 
expansion. Likewise, we have 

{A,{QA,{Q) ^ {a['\QA^^\Q) = 4(Ci(te)C2(te)) ^ 4(Ci(t)C2(t)) (10.6) 

where the last (approximate) equality is justified as long as the perturbative expansion (in 
fields) holds, since the difference between t and i is of first order in the fields. 
We can do the same for the 3-pt correlators: 

(A^) ~ (A^^U^^U^i)) + 3(^(1)^(1^(2)) _ g^^3^ ^ {Cdndn) . (10.7) 
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Notice that {('^dndn) ~ {dndny ~ 0(e while (C^) ~ 0{e thus if we restrict ourselves 
to the leading order in slow-roll, this term can be safely neglected. It follows that 

{Ai{ie)A2{ie)As{ie)) ~ 8 (Cl (te)C2 (te)C3 (^e) ) (10.8) 

^ 8(Ci(t)C2(t)C3W) + ^{Ut)^LAt)C2{t)Ut)) + perms. 

^ 8(Ci(t)C2(t)C3(t)) + 6 ■ -^m) Ci(t)C2(t)C3(t)) + perms. , 

where we have used that in the long wavelength limit ( ~ ('/a ~ ^c^eH (. In the last line, 
the second term and its permutations are negligible at the leading order in slow roll, since 
eiC) ~ e{Cr ~ 0(6-1) ^hile (C^) ~ 0{e-'). 

Hence, as long as we focus only on the leading contribution in slow roll, the first issue 
mentioned above can be easily resolved: the 2-pt and 3-pt correlators of scalar perturbations 
in UG are related to those of (2 times) ( in SFSG. Not surprisingly, similar relations for 
tensor perturbations hold if we apply the same logic. 

As for the second issue. It can be shown [281 129] that during the post inflationary 
era, when the matter content of the universe is in the form of a perfect fluid, the scalar 
perturbation A is adiabatic in the long wavelength limit, i.e. it is a constant at nonlinear 
level as long as it stays outside the horizon. However, unlike other inflationary models 
where there exists a conserved scalar mode in the long wavelength limit during inflation, 
the scalar perturbation A in our model evolves slowly outside the horizon, in the sense that 
A ~ A^^^ = 2( oc (— r)3^T^. Therefore, after the rapid transition from solid phase to perfect 
fluid phase, rather than staying at its initial value, A{te), it approaches its eventual constant 
value. However the relative difference between the two is only of order e: the slow time- 
dependence of A during inflation means that right after reheating the initial condition for the 
velocity is roughly A(te) ~ e ■ HA{te)- Since then, A decreases like 1/a^, thus making A{t) 
approach its asymptotic value in a few Hubble times, during which A{t) moves by ~ eA{te). 
At the leading order in slow-roll, we can neglect this difference. Notice that this effect cannot 
change the tilts that we have computed: all modes of interest are outside the horizon during 
reheating and during the phase when A relaxes to its asymptotic value. As a result, this 
small correction of order e to the value of A is the same for all modes, i.e., independent of 
k. By applying the same logic, we can reach the same conclusion for the transverse traceless 
tensor perturbation Dij, which is not conserved in the long wavelength limit during inflation, 
but approaches an asymptotic value in the post-inflationary fluid phase in a similar manner 
as its scalar counterpart. 

In conclusion, at the order we are working, we can take the correlation functions for ( and 
7 that we have computed during inflation in SFSG, evaluate them right before reheating, 
and obtain in this way good approximations to the corresponding correlation functions in 
UG in the post-inflationary phase. In particular, even though our scalar perturbations are 
not adiabatic during inflation, at reheating they get converted to adiabatic ones, with the 
same asymptotic constant value of ( (up to 0(e) corrections) as they had at reheating. 
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11 Summary and outlook 



Our model differs drastically from more standard ones in its symmetry breaking pattern. 
In particular, time-translations are not broken: there are physical "clocks" — i.e., time- 
dependent gauge-invariant observables — , but they inherit their time-evolution from the 
metric, not from the matter fields. As a result, the systematics of the EFT for the associ- 
ated Goldstone excitations is completely different than the standard effective field theory of 
inflation. This has far reaching implications, some of which are directly observable. 

The observational predictions of our model can be summarized as follows: (i) ^ nearly 
scale-invariant spectrum of adiabatic scalar perturbations, in agreement with observations. 
(a) A nearly scale-invariant spectrum of tensor modes, with a slight blue tilt; the tensor- 
to-scalar ratio r ~ ec\ ranges from somewhat smaller than in standard slow-roll inflation, 
for ultra-relativistic longitudinal phonons with c\ ~ 1/3, to tiny, if they are non-relativistic. 
(in) A scalar three-point function with a novel shape — peaked in the squeezed limit, with 
non-trivial angular dependence on how the limit is approached — and a potentially very large 



1 1 



amplitude, as big as /nl ~ -t 



A fourth prediction, whose full analysis we leave for future work, is the presence of 
vector modes. Unlike usual inflationary systems, our solid features transverse phonons, 
which get excited during inflation. In Appendix [C] we show that the vector-to-scalar ratio 
at the end of inflation scales as (c|_/c|.)^/^, which, depending on the value of c^,, ranges from 
negligible to roughly 6%. After reheating, the universe is dominated by a perfect fluid, and 
these vector modes should decay in the usual fashion, thus leaving no detectable imprint on 
the CMB. However, they interact at non-linear order with scalar and tensor perturbations 
during inflation and at reheating, thus affecting in principle their (higher point) correlation 
functions. 

Also for future work we leave a more thorough understanding of the time-evolution of 
our perturbations. One should be able to get model-independent information by running a 
separate universe-type argument involving anisotropic background solutions. In particular, 
this approach might elucidate whether and how super-horizon perturbations get affected by 
local processes at reheating. We showed in sect. 10 that a quick phase-transition triggered 
by det B^^ conserves C, and the tensor modes. Moreover, we showed that for the two- and 
three-point function, there is no difference between using unitary gauge time and SFSG time, 
which implies that we would get the same post- inflationary correlation functions if reheating 
were in fact triggered by another "clock", say [B]. This gives us confidence that, at least for 
these correlation functions, our predictions are robust. However it leaves open the possibility 
that we would get different predictions if the reheating phase lasted longer, for an Hubble 
time or more. 

There is also a number of generalizations of our model that we feel deserve being studied. 
The first, would be to promote our sohd to a 'super-solid'. In physical terms, a super-solid 
is a solid harboring a superfluid |30]. In our field-theory, symmetry-breaking terms, it is a 
system of four derivatively coupled scalars, 0° and 0^, with a shift-symmetry on 0°, our solid 
symmetries on the 0^'s, and a state that spontaneously breaks all space-time translations as 
follows [ID]: 

(00)= /it, (0^)=ax^ (11.1) 
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This system then combines the symmetry-breaking pattern of standard inflationary models 
with ours, and provides the minimal ingredients — the analog of 'single- field' — for studying 
the consequences of breaking all space-time translations during inflation, in a way that is 
consistent with a residual physical homogeneity and isotropy. Cosmological systems of this 
sort have been considered briefly in [8]. Notice that such a model will have two scalar 
perturbations — roughly speaking, excitations of 0" and longitudinal excitations of 0^, al- 
though in general they will be mixed, like it happens for instance in finite-temperature su- 
perfluids [21], where the two resulting modes go under the names of first and second sound. 
Notice also that, for inflationary purposes, one does not need to assume that the shift- 
symmetry on (j)^ be exact. One can assume an approximate shift-symmetry, with symmetry- 
breaking couplings suppressed by slow-roll parameters — in which case, to be consistent with 
the literature, one should refrain from calling such a system a super-solid. 

Another generalization that one should consider, is demoting our solid from isotropic 
to crystalline. We say "demoting", because it would entail lowering the degree of internal 
symmetries acting on our 0^ fields, from the full 50(3) group to one of its discrete subgroups, 
e.g. the cubic symmetry group. Then, in the background configuration with (0^) = ax^ , only 
that particular discrete subgroup will be preserved (now as a linear combination of internal 
and spatial rotations), and one could have interesting, observable deviations from isotropy. 
One of course needs to make sure that these deviations are absent from the background 
evolution and from the spectrum of perturbations — which we know to be isotropic with 
very good accuracy — but this might be automatic for certain subgroups of S0{3). As an 
example, consider the background evolution. The potential anisotropy in it is determined 
by the tensor structure of Tij. Let's assume for definiteness that our discrete subgroup is 
the cubic group — by which we also mean the individual inversions along the sides — and let's 
align the sides of the cube with x, y, and z. Then, Tij has to be invariant under permutations 
of X, y, and z, and under the individual parities x — ?■ —x, etc. The only two-index tensor 
with these properties is 



which is nothing but 5^^ . That is, at the two-index level, cubic symmetry accidentally implies 
full 5*0(3) invariance. For the spectrum of perturbations to be accidentally isotropic in the 
same fashion, one needs the phonons' quadratic action to be accidentally isotropic. This 
involves now four-index background tensors, contracted with the phonons' derivatives. So, 
the general mathematical question is, what are the discrete subgroups of 50(3) whose invari- 
ant two-index and four-index tensors are all isotropic? If one finds one or more subgroups 
with this property, and if their six-index invariant tensors are not all isotropic, one has a 
model of inflation, that, because of symmetry, has an isotropic background and spectrum of 
perturbations, and a (potentially) maximally anisotropic three-point function. 

Finally, it would be interesting to run a dedicated numerical analysis of CMB data for 
our specific three-point function template, given its small overlaps with the more standard 
templates that have been considered so far. 

We hope to address all these questions in the near future. 



X 
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Appendix 

A Time dependence of background quantities 

In order to solve the classical equations of motion for scalar and tensor perturbations, we 
need know the explicit time dependence of quantities such a(r), H{t), e(r), . . . ; the goal of 
this section is obtaining this time dependence. For the computations we are interested in, 
it suffices to derive these temporal functions up to the first order in slow roll. To make the 
notation lighter, we will mostly drop the r argument: a(r) — > a, etc. Primes will denote 
derivatives with respect to r. 



Recall the definition of the first slow roll parameter e, (3.6), and rewrite it as 

H' d f I 



aW^ dr \ aH 



(A.i: 



Integrating the above equation once and choosing some suitable additive constant one 
has 

^ = -(l-e,)r + 0(e^) (A.2) 

where the subscript "c" denotes evaluation at some reference conformal time Tc, which we find 
most convenient to choose to be the (conformal) time when the longest mode of observational 
relevance today exists the horizon, i.e. |ci^c^minT"c| — |cL,cTc-f^today| = 1- 

The reason e(r) is being treated as a constant in integration is that the higher order 
terms in the Taylor series of e(r) = e(rc) + e'(rc)(r — Tc) + . . . are suppressed by powers of 
slow roll, for instance e'{T^Tc ~ O(e^). Of course this also depends on the choice of reference 
time; we don't want the perturbative expansion in slow roll to be contaminated by large 
values of t/tc — 1. Since |r(t)| is a decreasing function during infiation, and +oo > —r > 0, 
we ought to choose early times (like Tc) as the reference. 

Using the definition of the Hubble parameter, we can extract the time dependence of the 



scale factor a(r) from the above equation (A.2): 



aH=- = — =^ a(T) = aJ-] + O .(e^) (A.3) 

a T \Tr 



'^The integration constant is chosen by demanding a(T) ^ a(rc), for t/tc — )• 0. 
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Furthermore, we obtain 



H(t) 



acTc 



(A.4) 



Finally, the time dependence of e, cl and ct can be revealed by invoking the definitions of 
other slow roll parameters. For instance, 



^- = aHr] = -"^ + 0(e' 
e r 



e r 



Similarly we obtain 



CL r 



ct,< 



(A.5) 



(A.6) 



Notice that, because of the all-order relation between and c| of footnote 4, ct,c and Mc 
are not independent parameters — they can be expressed in terms of cl,c and of the slow-roll 
parameters. The equations (A.3)-(A.6) are frequently used in solving the classical equations 
of motion. 



B Unitary gauge vs. spatially flat slicing gauge 

In performing calculations throughout this paper two gauge choices are particularly useful: 

• Spatially Flat Slicing Gauge (SFSG) is defined by setting to zero the scalar and vector 
perturbations in Qij, i.e. by imposing 

gij = a{tfexpjij , (B.l) 

where 'jij denotes the transverse traceless tensor mode, satisfying 

7n = , dm, = . (B.2) 

Then the fiuctuations in our 0^ scalars are unconstrained: 

0^ = + TT^ (B.3) 

The three tt^x) fields can be split into a transverse vector and longitudinal scalar as in 
(2.15), according to their transformation properties under the residual rotation group. 
This gauge choice is of particular convenience for computations of the two- and three- 
point functions because in the demixing (with gravity) limit the vr Lagrangian will 
contain all the scalar (or longitudinal) and vector (or transverse) degrees of freedom. 

• Unitary Gauge (UG) is defined by setting to zero the fiuctuations in the 0^ fields and 
in the "clock" field: 

(f)^ = x\ det{B^^) = a(t)-^ . (B.4) 
Then the spatial metric is unconstrained. And can be parameterized in general as 

g.. = a(tf exp(v4 5ij + didjX + diCj + djCi + Aj) , (B.5) 
where diCi = and Dij is transverse-traceless. 
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From the above form of the metric it seems that in UG there are too many degrees of 
freedom; there is an extra scalar in addition to the scalar, transverse vector, and transverse 
traceless tensor that we expect. However, when the metric is expressed in terms of the ADM 



parameters defined by (5.1), the second condition in (B.4) can be rewritten as 

+ V^x = log(l - N'Ni/N'^) . (B.6) 
As A^' and A^ can be expressed in terms of A, x, Cj, and Dij by solving the constraint 



equations given by (5.4), (5.5) we can see that (B.6) implies that the two scalar functions 
A{x) and x{^) ci-re not independent in UG. Hence the dynamical d.o.f. in question can be 
chosen to be A{x) (the scalar mode), Ci{x) (the transverse vector mode), and Dij (the 
transverse traceless tensor mode). The number of which matches our physical intuition and 
properly agrees with SFSG. UG is particularly useful in following our degrees of freedom 
through the reheating surface. 

As we find it convenient to calculate correlation functions in SFSG, and yet utilize UG to 
most easily describe the surface of sudden reheating, we need to develop the transformation 
rules to go from one gauge to the other. Let's denote by {x^} the coordinate system in 
SFSG and by {x^} that in UG. A gauge transformation relating SFSG to UG is given by 

= + ^^{x), where 

e°(x) = -^9.7r^(x) + 0(7r2), e{x)=7i\x) + 0{7,') (B.7) 
and the scalar perturbations are related by 

A = ^9,7r* + 0{ti^) . (B.8) 

C Vector perturbations 

In this section we will derive the spectrum for the vector modes vr^ (in SFSG), in the same 
manner as for scalar perturbations and tensor perturbations. Let's begin by writing the 
vector modes as 

4(fc,t) = 5^el(^)7rT,,(fc,t), (C.l) 

A=± 

where the polarization vectors satisfy the transverse condition ki€^{k) = and they form an 
orthonormal and complete set e\{k)e\,{k)* = 6\y. As before, writing 



(fc, t) = 7r^'(A?, t) dx{k) + 7r^'(fc, t)* d\{~k) (C.2) 



the creation and annihilation operators obey the usual commutation relation [dx{k), dy{k')'^] = 
{2'KY5\\i5^{k — k') and the classical equations of motions for tt^' follow from varying the 



quadratic ttt action (6.3), which is given by 

+ ?,HN^ + AH'^e cItt!=^ = (C.3) 
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where the second one follows from (5.13) (after is expressed in terms creation and an- 
nihilation operators in the same manner as vr^). Eliminating tt^', we reach a second order 
differential equation for using conformal time and keeping terms only up to the first 
order in slow roll, it reads 



dr 



n: 



d 



(C.5) 



3c? 



where u = ct/Hct = + O(e^). The general solution for the equation above takes the 

form 

1— u, 



d 



SHil} {l + U,)cT,ck\T, 



(C.6) 



where /3 = — | — 2ec + ^r/c + Uc and z^y = f + l^c — c| + ^rjc- The constants S and J-" are 
determined via the initial condition, which is specified by 



lim N^{k,T) = lim 



AH^ae d 



rT^oo A;2 ^^"^T ' ^p^^ ^3/2' 

where in the second equality we have used the fact that the mode function of the canonically 
normalized transverse vector field vr^^JJ'/ should match that in the fiat space vacuum at early 
times, that is 



. y/AecrH 1 



-icj-(T) fcr(l+Uc) 



(C.7) 



_^ + iJmf) 
Therefore, J-" = and 



-1/2 



1 



-iCT{T)kT(l+Uc) 



8 



Ur 



T,d 



, ,ef(.v+i/2) 



2e MpiHa^ 



-2ec+'r}c/2+Uc 



+ 0{e 



(C.8) 



(C.9) 



Then, by (C.3), we can derive the expression for vr^'; in particular, its asymptotic behavior 
at very late time when the mode exits the horizon is given by 



lim TT^{k,T) 

fcT-s>0- 



3 T,c<^c 



-Ct,c k TcYl--' 



Vc 5uc 
2 2 



-3Hr 



^cMpic^lk^/"^ 



(C.IO) 

the time dependence of which is the same as the tensor (6.17) and scalar (6.30) perturbations. 
The vector-to-scalar ratio simply 

|2 ^^2 ^5/2 

+ 0{e). {C.U] 



lim 

fcr-S>0- 



Given the relation (3.18), this is always smaller than (1/3)^^^ ~ 6%, which is reached only 
in the ultra-relativistic case with c| ~ 1/3 and ~ 1, or equivalently when |Fy + Fz\ ^ 
0{e^)\F\. 
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D Full trilinear action in SFSG 

Expanding the Lagrangian to third order in fluctuations about the FRW background in 
SFSG we have after a straightforward but lengthy computation 



22 



(1) (1) (1) (1) (1) 



(1) . (1) (1) (1) (1) (1) 

(1) . / \ 
- 6N Ml.Hi - a\r- N'f - 4[n^n] + (1 - c^p^] - (1 + ^ - 2cl)[Il]A 

+ [lif {\Fxxxa-^ - §^{Fxz + FxY)a-^ + i^^z + f^Fy) 

+ [n][n2] + FxY)a~' - iFz - ^^Fy) 

+ [n][n^n] (2Fxxa-' + I{Fxz + FxY)a-^ - §Fz - |Fy) 

(1) (1) 

+ §^Fzm + (iFz + iFy) [n^nn] - l{Fy + Fz)a\r- iV')9,7r^(7T^- N^) 

+ a2[n](7T^- N^f {-2Fxxa~' + ^,{Fy + Fz)) | , (D.l) 

where 11 denotes the 3x3 matrix diix^ and [■ ■ ■ ] indicates the trace; for instance, [n-^IUI] = 

Now, as discussed in Sectionjsjone only needs to solve the constraint equations 5S/ 5N = 
and SS/SN^ = to linear order in perturbations. We therefore don't need to worry about 
terms in A^, A^* that are quadratic in the fluctuations contributing to the cubic Lagrangian. 



The solutions to these equations are given by (5.11), (5.12), and (5.13). 



In particular, we are interested in two separate limits for computing the three-point 
function. The first is the de-mixing regime where k ^ aHe^^^. In this limit, to lowest order 
in slow roll, we can effectively set 6N and Nl to zero. Furthermore, note that all terms 
that are not of the form 11^ (like the final [n]7r^ term) are explicitly suppressed by slow roll 
parameters. We are left with 

£3 = aitfM^^H^^\^^{dnf - Idndj-K^du-K^ - pTidj-K^d^Ti^ + Idju'd^Ti^dkii'^ , (D.2) 



where we neglected boundary terms. We have freely used the Friedmann equations (3.7), 
various definitions of slow roll parameters, and the the total derivative 

det(97r) = ie,,,Q„„5V'9%™9'=7r" = \ {[d-nf - ^[d-nWdi:^] + 2[d7:^]) (D.3) 

(which is a total derivative because of the e-tensor structure). 

^^Since we are after the three-point function for scalar perturbations, we ignore the interaction between 
the tensor mode 7 and the tt fields. 
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The second limit is in the strong mixing (with gravity) regime. This occurs when k <^ 
aHe^f^ . In this hmit we can write to lowest order in slow roll 



5N 


k d 


f- 


a dr 




Nl 


= TTl , 






= T^T ) 





ke ttl 



(D.4) 

(D.5) 
(D.6) 



where we have estimated the time dependence of vr via the explicit classical solution to the 



first order equation of motion given by (7.8a). When we insert the above expressions into 



the full cubic Lagrangian we see that all the terms involving these auxiliary fields are going 
to vanish, as the reccurring combination (tt* — iV*) vanishes to lowest order in slow roll and 
6N is explicitly of order e^dn). So, surprisingly, we see that to lowest order in slow roll 



we recover the same expression (D.2) for the cubic Lagrangian in the strong mixing limit. 
This is a convenient fact, as it allows us to effectively use the same expression for the cubic 
interactions during the whole inflationary phase in our calculation of the three-point function 
in Section [7l 



E Three point function integral for general triangles 

For momenta configurations that are not roughly equilateral, we need to be a little more 
careful in computing the integral (7.6), as there are many different regions in the integral, 
each of which will necessitate the use of different expressions of vrx(r, k) given by (7.8a) and 
dfSbl ). 

Assume for the general case that ki > k2 > k^ and define three separate times ri, T2, and 
^3 by 

- cikiTi = e , -CLk2T2 = e , -cik^T^ = e , (E.l) 

and for the observation time r require |r| < |ri| — that is, all modes are very long at the time 
of observation. 

Splitting the object in question into regions defined by the times given above we can 
write schematically: 

3 / „T-3 i'T2 /-ri l-T 

/(r;-oo) = YlTTL{T,ki) ■ I dr' + dr' + dr' + dr' 



X 



/ j=i 



l[nUT',k,) + C.C. (E.2) 



We can do each integral separately, and will find that provided that the the triangle does 
not become "too squeezed" (which will become precise in a moment) only the first integral 
contributes to leading order in slow roll. For each integral there is, of course, going to be an 
overall H^^ULi |i3fcJ'(-CLA;.r)4^+^ as |r| < 
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E.l Integrating from — oo to r^: 

Ignoring for a moment the overall factor, we need to compute 



dr' 



T 



/4 



n 1 



i=l 



k.r' - Iclkfr'^ \ e*^^^'^' + c.c = h • 



(E.3) 



In order to ensure convergence of the integral and project onto the right vacuum the integral 
is computed over a slightly tilted contour, that is r' — )■ (1 — ie){T' + rs), with e > 0, and the 
limits of integration are from — oo to 0. Performing the integral and expanding in we have 
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3780 



{-3{kl + kl + kl) + l{klkl + 5 perms)) ■ r| + 0{tI) , (E.4) 



where U{ki^ /c2, k^) is the scale invariant function given by (7.13). Now, as = —e/cLk^ we 



can see that provided ki/k^ < {y/e) ^ only the zeroth order term in contributes to leading 



order in slow roll. That is, for a not "too squeezed" triangle the first integral in (E.2) is 
given by 

clkik2k3U{ki,k2,k3) ^ ,2 , -a a ^ ^ /-p 

~ 97ff2 11 l^fe^r X (1 + C(e)) , provided ^ < • (E-5) 



i=l 



E.2 Integrating from to 

Again, dropping the overall factor, first note that 

2 



T"3 



dr' 



-ciksTs) 



n/4 



-CLk^r' 



i=l 



nil- tCLhr' - '-clk^r'^ ) e^'^-'^^' + c.c 



This final form of the integral is straightforward to compute and yields 

Cl{kl + kl){k2-k3){k2 + k3) 



C.C 



(E.6) 



Abklkl 



(E.7) 



which, provided that /ci/fcs < (Ve) ^ is first order in slow roll. That is, for a not "too 



squeezed" triangle the second integral in (E.2) is given by 



^n|S.J^xO(e), provided ^<^. 



(E.8) 
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E.3 Integrating from r2 to ri: 

The exact same logic follows for this piece. The meat of the integral, 



£ dr (l - iCLk.r' - \clklT'^^ e^-^^^' + c.c ^ h , (E.9) 



yields 

^3 / 1.2 



45 M k 



I, = -!^ki[l--^\-e' + 0{e^). (E.IO) 



Which, provided that ki/k2 < (a/c) ^ is first order in slow roll. And so, once again, for a 
not "too squeezed" triangle the third integral in (E.2) is given by 



n ^ ^(^) ' provided ^ < ^ . (E.ll) 



2=1 



E.4 Integrating from ri to r: 

As done in the main body of the text, we once can immediately see that this integral is going 
to be of the form: 



i=l 



dr' {-i)g{T') + c.c. 



(E.12) 



where gir) is pure real. And thus, the final integral in (E.2) simply vanishes. 



E.5 Summary: 

And so, for not "too squeezed" triangles, /(r; —oc) is completely dominated by the portion 
of the integral where all the modes are in the de-mixed regime, that is: 

i 
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